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4. Prüfverteilungen 


4.01 a) 1,725 b) 2,896 

4.02 a) 2,086 b) 2,528 

4.03 Needs["Statistics‘Master«"] 

vert = StudentTDistribution[10]; 

Plot[PDF[vert,x],{x,-5,5>,GridLines->Automatic,Frame->True]; 



4.04 a) 31,41 b) 20,09 c) 43,82 
4.05 a) 8,26 b) 2,18 c) 17,19 
4.06 a) 20,48 b) 66,77 c) 39,38 
4.07 vert = ChiSquareDistribution[4J; 

Plot[PDF[vert,x],{x,0,10},GridLines->Automatic,Frame->True ]; 



4.08 a) 9,36 
4.09 a) 2,20 
4.10 a) 0,361 


b) 7,39 
b) 6,62 
b) 0,145 


c) 5,47 
c) 2,96 
c) 0,0775 




4. Prüfverteilungen 


4.11 vert = FRatioDistributionf 6,10]; 

Plot[PDF[ver t,x],{x,0,10},GridLines->Automatic,Frame->True]; 



4.12 Wir laden die Hilfsdatei PROBABIL.MTH 

IF(t < 0, 0.5 (1 - STUDENT(t, 10)), 0.5 (1 + STUDENT(t, 10))) 
IF(x < 0, 0, CHI_SQUARE(X, 4)) 

IF(x < 0, 0, 1 - F_DISTRIBUTION(x, 6 , 10)) 
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5. Zufallsstreubereiche und Vertrauensbereiche 
bei diskreten Verteilungen 

5.01 (1) a) 1 s X s 9 b) 0 * x s 7 

c) 1 s x s 8 

(2) a) 9 b) 7 

C) 7 

(3) a) 2 b) 1 

c) 1 

5.02 Larson-Nomogramm: Seite 9 

5.03 a) 0 s x s 7 b) 2 s x s 11 c) 6 * x s 9 

5.04 a) 6 b) 10 c) 18 

5.05 a) 1 b) 2 c) 7 

5.06 Thorndike-Nomogramm: Seite 10 
5.07 Ergebnisse: 

(1) a) 0,77 % s p s 20,8 % b) 1,38 % s p < 12,3 % 

c) 2,3 % s p s 11,2 % 

(2) a) 18,4 % b) 11,1 % c) 10,3 % 

(3) a) 1,1 % b) 1,7 % C) 2,7 % 

5.08 Bi: n = 100 x = 5 a = 0,05 1,6 % s p < 11,3 % 

x F 

5,09 p unten = n-x+l+xF u mit F u " F 2x, 2(n-x+l); a 
(x + 1) F q 

p oben “ n - x + (x + 1) F q mit F o “ F 2(x+1), 2(n-x); 1-a 

5.10 (1) a) 0,24 s 4 s 7,22 b) 1,09 s 4 s 10,24 

C) 3,45 s 4 ^ 15,76 

(2) a) 6,30 b) 9,15 c) 14,43 

(3) a) 0,36 b) 1,37 c) 3,98 

5.11 Thorndike-Nomogramm: Seite 11 

5.12 Po: x = 3 a = 0,05 0,62 s 4 s 8,77 


5.13 Po: n 


= 0,01 


0 s 4 s 1,06 



5. Zufallsstreubereiche und Vertrauensbereiche 


5.02 Daten von 5.01 c) 

























































































6. Statistische Tests 


6.01 Webstuhl A: 5 Ballen mit insgesamt 10 Fehlern 
Webstuhl B: 5 Ballen mit insgesamt 4 Fehlern 

1. Schritt: Die Nullhypothese lautet: m a = 

2. Schritt: a = 1 % 

3. Schritt: Berechnung der Prüfgröße 


4. Schritt: Überprüfung der Voraussetzung 
Beide "zu erwartenden" Werte (x a ., = : 


Daher ist die Voraussetzung für den x -Test erfüllt. 


6. Schritt: Die Prüfgröße ist kleiner als der Tabellenwert. 

7. Schritt: Interpretation 

Obwohl auf den 5 Ballen des Webstuhls A mehr als doppelt 
so viele Fehler wie auf den 5 Ballen des Webstuhls B 
gefunden wurden, besteht kein signifikanter Unterschied 
zwischen den beiden Webstühlen hinsichtlich der mitt¬ 
leren Anzahl von Fehlern. 

6.02 1. Charge: 16 Risse in der Prüfeinheit 

2. Charge: 9 Risse in der Prüfeinheit 

3. Charge: 21 Risse in der Prüfeinheit 


Nullhypothese: = ß 2 = ü 3 


ü = = 15,33 


1. Charge 

2. Charge 

3. Charge 

*prüf = ^ J 


beobachtet - erwartet) _ 
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Alle 3 "zu erwartenden" Werte x x th = x 2 th : 
sind a 5. Daher ist der z 2 -Test zulässig. 

1 o,95 — ^>991 < X+ 


^tab ‘ 


*prüf 


Interpretation: Die Abweichung der einzelnen Stichproben aus 
den drei Chargen untereinander kann auch zufällig sein. 

Ein Unterschied zwischen den drei Chargen hinsichtlich der 
mittleren Anzahl der Risse kann statistisch nicht nachge¬ 
wiesen werden. 


6.03 2. Charge: 9 Risse in der Prüfeinheit 
3. Charge: 21 Risse in der Prüfeinheit 

Nullhypothese: H Q : = ß 2 = ß a = 5 % £ = 15 



beobachtete AvF 

erwartete AvF 

2. Charge 

9 

15 

3. Charge 

21 

15 


'‘prüf 

z 2 = 
*tab 


Die Voraussetzung für den z -Test ist erfüllt. 
= 3,841 


*prüf 


*tab 


Interpretation: Der Unterschied zwischen 2. und 3. Charge ist 
signifikant. 

Beachten Sie: 

Beim Vergleich von allen 3 Chargen war der Unterschied nicht 
signifikant. 

Der Unterschied zwischen 2. und 3. Charge (ohne den "ausglei¬ 
chenden Einfluß" der 1. Charge) war signifikant. 

Je nach Fragestellung erhält man eine andere Antwort. 


6.04 Nullhypothese 

a) VB(m a ): 

U B * VB ^ a > 

b) VB(/i a ): 

U ß e VB(m a ) 

c) VB(m a ): 

M ß (S VB(4 a ) 


48,88 s 4 a s 49,32 
=» Hg wird verworfen 
12,41 s (i A s 12,45 
=» H q wird nicht verworfen 
182,67 s u A s 183,33 
=» Hg wird verworfen 
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6.05 


Nullhypothese H Q : m a = 

Wir müssen mit der t-Verteilung arbeiten. 
Es gilt: M A = x A ? t n _ 1;1 _ a/2 • s / Vn 

a) VB(H A ): 48,91 s s 49,29 

/i ß * VB(m a ) =* H Q wird verworfen 

b) VB(4 a ): 12,42 £ U A * 12,44 

M ß e VB(m a ) =» H Q wird nicht verworfen 

c) VB(m a ): 182,71 s m a s 183,28 

Mg «s VB(m a ) =» Hg wird verworfen 


6.06 Grundgesamtheit A: <r = 0,2 g, daraus Stichprobe mit 

= 499,86 g 

Grundgesamtheit B: m = 500 g 

Zweiseitige Fragestellung Nullhypothese H Q : m a = M ß 

Wir wählen a = 5 % , 1 % und 0,1 % 

Zweiseitig abgegrenzter (1 - a)-VB bei bekanntem er: 


1 - a 

u l-a/2 

"un 

4 ob 


95 % 

1,96 

499,77 

499,95 


99 % 

99,9 % 

2,58 

3,29 

499,74 

499,71 

499,98 

500,01 

M ß = 500 


6.07 
6.08 
6.09 
6.10 


Mg = 500 g liegt weder im 95 %-VB noch im 99%-VB von m a - 
Es liegt lediglich im 99,9 % von m a * 

Interpretation: 

Die Nullhypothese H Q : m a = M ß wird bei a = 5% und bei 
a = 1% verworfen. 

Eine Abweichung vom Sollwert ist signifikant, d.h.stati¬ 
stisch nachgewiesen. 

Die Maschine ist - sehr wahrscheinlich - nicht richtig 
eingestellt. 


7 

4,26 

H o 

wird 

verworfen 

6,8 

3,1 

H 0 

wird 

verworfen 

0,496 

3,24 

H 0 

wird 

nicht verworfen 

1,21 

2,76 

H 0 

wird 

nicht verworfen 




6. Statistische Tests 


6.11 0,84 4,26 H q wird nicht verworfen 

6.12 •) *p rüf - 21.14 * 31 ; 0.95 * 45 

H q : <t a = o-g wird mit a = 0,05 nicht verworfen 

*» *p rUt - 68,4 231,0,95 ‘ 66 ’ 3 

H q : ct a = o-g wird mit a = 0,05 verworfen 

6.13 a) 


xquer= 82,166667 


Max= 

Mln= 

AnzahU 


89 

76 


2,6691373 
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6 . Statistische Tests 


Beispiel von Seite 76 mit MATHEMATICA. 

Needs["Statistics'Master'"] 

(* Chi-Quadrat-Anpassungstest für Normalverteilung" *) 
nnv=NormalDistribution[0,1]; 

n » { 1, 2, 3, 5, 6, 10, 14, 

16, 19, 17, 13, 6, 7, 5, 1}; 

x = {14.8, 14.9, 15, 15.1, 15.2, 15.3, 15.4, 

15.5, 15.6, 15.7, 15.8, 15.9, 16, 16.1, 16.2>; 

cob = {14.75, 14.85, 14.95, 15.05, 15.15, 15.25, 15.35, 15.45, 

15.55, 15.65, 15.75, 15.85, 15.95, 16.05, 16.15, 16.25}; 

alpha = 0.05; 

k=Length[n] 15 

sn=Sum[n[[j]],{j,l,k}]; 

xq=n.x/sn 15.5656 

s=Sqrt[(n.(x A 2) - sn xq A 2)/(sn-l)] 0.288797 

u=Table[ (cob[ [j]]-xq)/s, {j,l,k+l}]; 
g=Table[CDF[nnv,u[[j]]],{j,l,k+l}]; 
dg=Table[g[[j]]-g[[j-l]], {j,2,k+l}]; 
er=Table [sn dg[[J]],{j,l,k>] 

{0.529764, 1.23891, 2.57294, 4.7452, 7.77174, 11.3038, 14.6007, 
16.7482, 17.0613, 15.4348, 12.4004, 8.84741, 5.60582, 3.15429, 
1.57616} 

Wir überprüfen die Voraussetzungen: 

Alle "zu erwartenden" Häufigkeiten, also alle Elemente des 
Vektors er, müssen mindestens 1 sein, höchstens ein Fünftel 
der Werte darf < 5 sein. 

Wir fassen die ersten 3 Klassen zu einer neuen Klasse mit 

= 14.9 und n^ = 6 zusammen, ebenso die beiden letzten Klassen 
zu einer neuen Klasse mit x 12 = 16.15 und n 12 = 6. 


Needs["Statistics‘Master *"] 

(* Chi-Quadrat-Anpassungstest für NormalVerteilung" *) 
n = { 6, 5, 6, 10, 14, 16, 

19, 17, 13, 6, 7, 6}; 

X = {14.9, 15.1, 15.2, 15.3, 15.4, 15.5, 

15.6, 15.7, 15.8, 15.9, 16, 16.15}; 

cob={14.75, 15.05, 15.15, 15.25, 15.35, 15.45, 

15.55, 15.65, 15.75, 15.85, 15.95, 16.05, 16.25}; 
alpha = 0.05; 

k=Length[n] 12 

sn=Sum[n[[j]],{j,l,k}]; 

xq=n.x/sn 15.5656 

s=Sqrt[(n.(x A 2) - sn xq A 2)/(sn-l)] 0.294945 

u=Table[(cob[[j]]-xq)/s, {j,l,k+l}]; 

g=Table[CDF[NormalDistribution[0,l],u[[j]]], {j, l,k+l}]; 
dg=Table[g[[j]]-g[[j-l]], {j,2,k+l}]; 
er=Table[sn dg[[j]],{j,l,k}] 

Die Voraussetzungen für die Anwendung des x 2 -Tests sind erfüllt. 
Wir dürfen weiterrechnen. 
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prüfvert=Sum[(n[[j]]-er[[j]] ) A 2 / er[[j]],{j,l,k}] 
tabellenwert= Quantile[ChiSquareDistribution[k-3],1-alpha] 

{4.6722, 4.89811, 7.86207, 11.2615, 14.3947, 16.4196, 16.7139, 
15.1825, 12.3073, 8.90296, 5.74718, 5.01263} 

2.96673 

16.9190 

x priif = 2,97 < x tab ~ 16,9 * H 0 wird nicht verworfen. 


Wir rechnen nun den zweiten Durchgang mit DERIVE. 

Die Berechnung von s ist heikel. Eine Summe von Quadraten von sehr 
kleinen Zahlen, und daraus die Quadratwurzel bringt Probleme, die 
wir bereits im 2. Jahrgang besprochen haben (Wahl eines vor¬ 
läufigen Mittelwerts, usw.). 

Bei der üblichen Genauigkeit erhalten wir s = 0,295072 
In der Hilfsdatei PROBABIL.MTH stellt uns DERIVE die Funktion 

CHI_SQUARE(X2, V) 

für die Berechnung der Verteilungsfunktion der ar-Verteilung 
zur Verfügung. 

18: PrecisionDigits := 12 
19: Notation := Decimal 
20: NotationDigits := 6 

21: n := [ 6, 5, 6, 10, 14, 16, 

19, 17, 13, 6, 7, 6] 

22: x := [14.9, 15.1, 15.2, 15.3, 15.4, 15.5, 

15.6, 15.7, 15.8, 15.9, 16, 16.15] 

23: cob := [14.75, 15.05, 15.15, 15.25, 15.35, 15.45, 

15.55, 15.65, 15.75, 15.85, 15.95, 16.05, 16.25] 

24: a := 0.05 

25: k := DIMENSION (x) 12 

27: sn := SUM(ELEMENT (n,j), j, 1, k) 125 

29: xq := (n . x)/sn 15.5656 

31: xx := VECTORJELEMENT (X,j) A 2, j, 1, k) 

32: [222.01, 228.01, 231.04, 234.09, 237.16, 240.25, 237.16, 
240.25, 243.36, 246.49, 249.64, 252.81, 256, 260.822] 

33: s := /( (n xx - sn xq A 2) / (sn - 1)) 0.294945 

35: u := VECT0R((ELEMENT(cob,j ) - xq) / s, j, 1, k+1) 

36: [-2.76526, -1.74812, -1.40907, -1.07002, -0.730983, -0.391937, 
-0.0528912, 0.286155, 0.625201, 0.964247, 1.30329, 1.64233, 
2.32043] 

37: g := VECT0R(NORMAL(ELEMENT(#36, j )) , j, 1, k+1) 

38: [0.00284386, 0.0402214, 0.0794063, 0.142302, 0.232394, 

0.347552, 0.478909, 0.612620, 0.734080, 0.832539, 0.903762, 
0.949740, 0.989841] 

39: dg := VECT0R( ELEMENT(#38, j) - ELEMENT(#38, j-1), j, 2, k+1) 
40: [0.0373775, 0.0391849, 0.0628965, 0.0900916, 0.115157, 
0.131357, 0.133711, 0.121460, 0.0984585, 0.0712236, 

0.0459774, 0.0401010] 

41: er := VECTOR(sn ELEMENT(#40, j), j, 1, k) 

42: [4.67219, 4.89811, 7.86207, 11.2614, 14.3947, 16.4196, 

16.7138, 15.1825, 12.3073, 8.90296, 5.74718, 5.01263] 

43: pw: = SUM((ELEMENT(n, j )-ELEMENT(er,j)) A 2/ELEMENT(er ,j ) , j, k) 
44: 2.96673 

45: chq := CHI SQUARE(2.96673, k - 3) 


0.0344006 
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Beachten Sie: 

DERIVE hat keine Funktion zur Berechnung des Quantils x^ a 

Wir müssen daher einen anderen Weg gehen. 

Wir berechnen G (Zp rüf>k _ 3 )• 

Wenn dieser Wert kleiner als 1 - a ist, dann wird unsere Null¬ 
hypothese H q : Es liegt eine NV vor, nicht verworfen. 

Wenn dieser Wert größer als 1 - a ist, dann wird unsere Null¬ 
hypothese H q : Es liegt eine NV vor, verworfen. 

0,0344006 <0,95 => H Q wird nicht verworfen. 

Beachten Sie: 

Für die meisten Tätigkeiten gibt es optimale Werkzeuge. 

Das Durchrechnen unseres Beispiels mit MATHEMATICA und DERIVE ist 
lediglich als Demonstration dieser Systeme gedacht. 

In der nächsten Aufgabe wird gezeigt, wie man Aufgaben dieser Art 
elegant mit EXCEL, einem Tabellenkalkulationsprogramm, löst. 
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6.14 a) 


Max= 89 

Mln= 76 

Anzahl= 60 


xquer= 82,16667 
8= 2,669137 
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Klassen 6 und 7 Zusammenlegen: 
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Summe= 7,9437 

Freiheitsgrad f = k - 3 = 3 

Chl-Quadrat-Prüfgröße = 7,9437 

Chl-Quadrat-Tabellenwert = 7,8147 

(zu 95 %) 


Nullhypothese verworfen, Daten sind nicht normalverteilt! 











22 


6. Statistische Tests 


6.14 b) 


Max= 304 xquer= 281,9867 

Mln= 253 s= 9,558711 

Anzahl= 150 
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Klassen 1 - 3 Zusammenlegen: 
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Summe= 14,6915 


Chi-Quadrat-Prüfgröße = 14,6915 (Freiheitsgrad f = k - 3 = 8) 

Chi-Quadrat-Tabellenwert = 15,5073 

(zu 95 %) 


Nullhypothese kann nicht verworfen werden, 
Daten scheinen normalverteilt! 









6 . Statistische Tests 


6.15 b) y = 0,390699 X + 1233,8967 

c) Korrelationskoeffizient: 0,27984 

d) Standardabweichung: 13,9480 

s g = 22,2759 
s b = 0,3718 

e) X = 0,20044 y + 29,669 

f) y(50) = 143,4 



6.16 b) y = 0,780992 x + 0,801653 

c) Korrelationskoeffizient: 0,82286 

d) Standardabweichung: 1,32658 

S a = 1,29046 
s. = 0,19068 
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6 . Statistische Tests 


Wir bringen noch ein BASIC-Programm für dieses Problem. 

' Lineare Regression 

DIM x(50), y(50), v(50) 

CLS 

READ n: PRINT "n = "; n 

sx = 0: sy = 0: sxx = 0: sxy = 0: syy = 0: vqu = 0 
FOR i = 1 TO n: READ x(i), y(i) 

sx = sx + x(i): sy = sy + y(i): sxx = sxx + x(i) A 2 
sxy = sxy + x(i) * y(i): syy = syy + y(i) A 2 
NEXT i 

xquer = sx / n: yquer = sy / n 
ssxx = sxx - n * xquer A 2 

ssxy = sxy - n * xquer * yquer 

ssyy = syy - n * yquer A 2 

b = ssxy / ssxx: a = yquer - b * xquer 

PRINT 

FOR i = 1 TO n 

v(i) = a + b * x(i) - y(i) 
vqu = vqu + v(i) A 2 
NEXT i: best = 1 - vqu / ssyy 
PRINT "sx,..."; sx, sy, sxx, sxy, syy 
PRINT "ssxx,ssxy,ssyy"; ssxx, ssxy, ssyy 
PRINT 

PRINT "xquer 
PRINT "yquer 

PRINT "Summe der Fehlerquadrate 
PRINT "Bestimmtheitsmaß B 
PRINT "Korrelationskoeffizient r 
s = SQR(vqu / (n - 2)) 
sb = s * SQR(n / (sxx * n - sx A 2)) 
sa = s * SQR(sxx / (sxx * n - sx A 2)) 

PRINT "Fehlervarianzen: s-quadrat 

PRINT " sb-quadrat 

PRINT " sa-quadrat 

PRINT 

PRINT "Fehlerstandardabweichungen s 
PRINT " sb 

PRINT " sa 

PRINT 

PRINT " Regressionsgerade: y = "; b; "* x"; 

IF a >= 0 THEN PRINT " + "; a ELSE PRINT " - "; ABS(a) 

PRINT : PRINT 
f = n - 2 

PRINT "(1 - alpha)-Vertrauensbereiche" 

INPUT "alpha in % "; alpha 

PRINT "t("; f; alpha / 100; "zweiseitig) "; 

INPUT t 

PRINT "chi-quadrat("; f; ";"; alpha / 200; ") "; 

INPUT chi2 

PRINT "chi-quadrat("; f; ";"; 1 - alpha / 200; ") "; 

INPUT chil 

bu=b-sb*t: bo=b+sb*t: au=a-sa*t: ao=a+sa*t 
cu = (n - 2) * s A 2 / chil: co=(n-2) *s A 2/ chi2 

PRINT alpha; "%- VB von beta = [ "; bu; bo; "]" 

PRINT alpha; "%- VB von alpha = [ "; au; ";"; ao; "]" 

PRINT alpha; "%- VB von sigma-quadrat = [ "; cu; ";"; co; "]" 

DATA 10, 5,7, 4,3, 9,8, 10,10, 9,7, 6,4, 7,6, 6,4, 3,4, 5,5 


xquer, 
yquer 
vqu 
best 

SQR(best) 



















GridLines->Automatic, FrameLabel->{"p in % 



0 2 4 6 8 

p in % 


1.75 

1.5 

?£ 1.25 

1 

□ 0.75 
0.5 
0.25 
0 



p in % 


DERIVE 


Wir- laden die Hilfsdatei PROBABIL. 
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7.02 b) Needs["Statistics'Master'"] 

tl=Table[{100 p f 100 CDF[BinomialDistribution[100,p],3]>, 

{p f 0, 0.1, 0.001}]; 

gl=ListPlot[tl, PlotJoined->True, Frame->True, 

GridLines->Automatic, FrameLabel->{"p in %","Pa in %"}]; 

t2=Table[{100 p,p 100 CDF[BinomialDistribution[100,p],3]}, 
{p, 0, 0.1, 0.001}]; 

g2=ListPlot[t2, PlotJoined->True, Frame->True, 

GridLines->Automatic, FrameLabel->{"p in %","D in %"}]; 

100 

^ 80 

•- 60 

(0 

Q- 40 
20 
0 



t- 


[iöc 

fco 

1 


x 











x 







— 


0 2 4 6 8 10 

p in % 



p in % 
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7. Anwendungen der statistischen Methoden im QM 


7.02 c) Needs["Statistics'Master'"] 

tl=Table[{100 p,100 CDF[BinomialDistribution[125,p],4]>, 

{p, 0, 0.1, 0.001}]; 

gl=ListPlot[tl, PlotJoined->True, Frame->True, 

GridLines->Automatic, FrameLabel->{"p in %","Pa in %"}]; 
t2=Table[{100 p,p 100 CDF[BinomialDistribution[125,p],4]}, 

{p, 0, 0.1, 0.001}]; 

g2=ListPlot[t2, PlotJoined->True, Frame->True, 

GridLines->Automatic, FrameLabel->{"p in %","D in %"}]; 

100 

SR 80 

60 
co 

0- 40 
20 
0 

0 2 4 6 8 10 

p in % 


_c 

Q 


0 2 4 6 8 10 

p in % 

7.03 Needs["Statistics'Master'"] 

tl=Table[{100 p,100 CDF[BinomialDistribution[50,p],0]}, 

{p, 0, 0.1, 0.001}); 

gl=ListPlot[tl, PlotJoined->True, Frame->True, 

GridLines->Automatic, FrameLabel->{"p in % ,, ,"Pa in %"}]; 

100 

^ 80 
-- 60 
co 

0- 40 
20 
0 

0 2 4 6 8 10 

p in % 
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t2=Table[{100 p, 100 CDF[BinomialDistribution[50,p],1]>, 

{p, 0, 0.1, 0.001}]; 

g2=ListPlot[t2, PlotJoined->True, Frame->True, 

GridLines->Automatic, FrameLabel->{"p in %","Pa in %"}]; 



p in % 


Needs [ '• Stat ist ics 'Master'" ] 

tl=Table[{100 p,100 CDF[BinomialDistribution[80,p],2]}, 

{p, 0, 0.1, 0.001}]; 

gl=ListPlot[tl, PlotJoined->True, Frame->True, 

GridLines->Automatic, FrameLabel->{"p in %","Pa in %"}]; 
t3=Table[{100 p,p 100 CDF[BinomialDistribution[50,p],2]}, 

{p, 0, 0.1, 0.001}]; 

g3=ListPlot[t2, PlotJoined->True, Frame->True, 

GridLines->Automatic, FrameLabel->{"p in %","Pa in %"}]; 

100 

£ 80 
■E 60 
o- 40 
20 
0 

Show[gl, g2, g3]; 





f~50 

” 2 | 


N 

v_ 
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7.04 Aus G(t) = 1 - e -t / T erhält man durch Umformen: 
t = -T ln R 

t = -20000 ln 0,9 = 2107,2 Nach 2107 Stunden 

7.05 Aus G(t) = 1 - e -Xt erhält man durch Umformen: 

X = -ln(l-G)/t 

X = -ln(l - 50/10000)/10 = 5,0125410 -4 h" 1 
T = 1/X = 1995 h 

X = 5,013•10~ 4 h' 1 T = 1995 h 


7.06 R(200) = e -200 / 1995 = 0 ,90461 

B(200) = B q • R(200) = 10000 0,90461 = 9046,1 


9046 Stück 


7.07 Aus 
t = 


R(t) = e~ X erhält man durch Umformen 
-ln 0,95/(4,825 10 -6 ) = 10630,7 


t = -ln R(t) / X 

10631 Stunden 


7.08 R(2 Jahre) = 0,85 
Aus R(t) = 


erhält man durch Umformen b 


b = 


ln(-ln 0,85) 
ln(2 / 54,5) 


1,796 


ln(-ln R) 
ln(t/T) 


7.09 R(20) = 0,99 
a) b = 1 

T = -t/ln R = -20/ln 0,99 = 1989,98 T = 1990 h 

R(100) = e ' 100/1990 = 0,95099 

B(100) = B q R( 100) = 50000,95099 = 4754,95 4755 Stück 


b) b = 
T = 



T = 


20 

1 . 8 /-' 

✓-ln 0,99 


257,50 


T = 257,5 h 


*( 100 ) - e - UOO/^.«) 1 ' 8 , „, 83352 

B(100) = B q R(100) = 5000-0,83352 = 4157,5 4158 Stück 



8. Aufgaben zu den statistischen Methoden 
des Qualitätsmanagements 


8.01 P(fehlerfrei) = 0,9 • 0,97 • 0,95 • 0,93 = 0,771 
P(fehlerhaft) = 1 - P(fehlerfrei) = 0,229 

8.02 Hy(N = 30; d = 7; n = 5) P(x = 0) = g(0) = 0,2361 

8.03 Bi(p = 0,07; n = 20) 

P(X a 1) = 1 - G(1) = 1 - 0,5869 = 0,4131 

8.°4x U n = 32 g x ob = 38 g u un = - 0,8 u^-1,6 

P = 0,2119 + 0,0548 = 0,2667 Ausschuß = 26,7% 

8.05 a) x Qb = 52 mm u Qb = 2/3 P(x > 52mm) = 0,2525 = 25,25% 

b > X un = 47 ' 3 “» u un " - °' 9 

P(x < 47,3mm) = 0,1841 = 18,41% 

8.06 a) x un = 1800 h u uR =3 P(x a 1800 h) = 0,0013 = 0,13% 

b) X Qb = 1300 h U Qb = -2 P(x s 1300 h) = 0,0228 = 2,28% 

c) x x = 1400 h x 2 = 1700 h u i = - 0,5 u 2 = 2 

P(X X s x s x 2 ) = 0,819 = 81,9 % 

8.07 ii = X + U*tr U Q 99g = 3,09 => ß = 92,472 

“eins < 0,1 % 92,47 V s ß s 97,53 V 

8.08 Einzelwerte: ß = 104 g <r = 5 g 

Mittelwerte: ß- = 104 g er- = 5/V2Ö g 

x un = 101 g P(x > 101 g) = 99,64 % 

8.09 a) 99 % ZSB(x| ß = 500 W; <r = IW) : 497,42 W s x s 502,58 W 

b ) u o j99 5 = 2,576 <T = (x - ß) /u = 0,7764 

Probe: ZSB(x|u = 500; <r = 0,776): 498,0001 s x s 501,999 
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95 CI s x s 105 £1 
4,05 n i S i 11,22 n 


6 

8.12 0 bis 6 

8.13 Normalverteilung: 511,71 g i (i s 512,29 g 

0,275 g s <r s 0,730 g 

8.14 Binomialverteilung, 95%-VB(p): 0,78% s p s 10,57% 

8.15 Poisson-Verteilung, 95%-VB(n): 0,69 s ß a 3,15 

8.16 ß = 600 N/mm 2 a = 10N/mm 2 n = 10 x = 591,7N/mm 2 

a) 99%-ZSB(x) un : x a 592,6 N/mm 2 

591,7 * ZSB => mit a = 1 % kann angenommen werden, daß 
die Festigkeit geringer geworden ist. 

b) 99 %-ZSB(s) ob : s =s 15,52 N/mm 2 

15,1 e ZSB =» gleiche Qualität kann bei a = 1 % nicht 
ausgeschlossen werden. 

8.17 99%-VB(h a ): U a Ä 2,511 kW X = 2,441 kW 

2,5 e VB =» Hg wird nicht verworfen s = 0,0785 kW 

Aufgrund der Stichprobe kann man nicht sagen, daß der Liefe¬ 
rant eine zu hohe mittlere Leistung angegeben hat. 

8.18 n = 25 x = 97,003 mm s = 0,009 mm 

95%-VB (ß): 96,9993 mm s ß < 97,0067 mm 

Die Toleranzmitte 97 mm liegt in diesem VB. 

Mit a = 5% kann gesagt werden, daß der Mittelwert der Produk¬ 
tion mit der Toleranzmitte übereinstimmt. 


8.10 N(100; 7,5; 99%-ZSB(x): 

99%-ZSB(s): 

8.11 Poisson-Verteilung mit ß = 3 
95%-einseitig abgegrenzter ZSB(x): 

x un,eins 1 x ob,eins 



34 


8 . Aufgaben zu den statistischen Methoden 


8.19 95%-ZSB(s ß ): s ß s 1,454g 

99%-ZSB(s ß ): s ß s 1,546g s A = 1,52g 

99,9%—ZSB(s ß ): s ß =£ 1,650g 

s A = 1,52g liegt im 99%-ZSB(s ß ), aber nicht im 95%-ZSB(s ß ). 
Mit 1% < a < 5% können wir sagen, daß das er der Abfüllanlage 
größer als 1,25 g ist. 

8.20 99%-ZSB(s ß ): 0,108 mm * s ß s 0,299 mm 

s A = 0,167 mm liegt im 99%-ZSB(s ß ). 

H q : <r A = <r ß wird nicht verworfen. 

Es kann nicht nachgewiesen werden, daß sich die Standardab¬ 
weichung der Anlage durch die Wartung verändert hat. 

8.21 ü = x = 121,8 mm a = sT/a 8 = 0,05 mm/0,965 = 0,0518 mm 

x : OEG = 112,85 mm UEG = 121,75 mm 

x : OWG = 121,84 mm UWG = 121,76 mm 

s : OEG = 0,088 mm UEG = 0,0195 mm 

s : OWG = 0,078 mm UWG = 0,026 mm 

8.22 Binomialverteilung mit n = 125; x = c = 2 

a) p = 0,04 

P a = G(c|p;n) FZT34 G = 0,1196 P a = 12% 

b) p = 0,01 

P a = G(2|0,01;125) = 0,8693 1 - P fl = 0,1307 

Rückweisewahrscheinlichkeit = 13% 

8.23 4 pr = 1,0 P a = G(3|1) = 0,98101 P fl = 98 % 




9. Determinantenrechnung 

9.01 a) | J ~g | = 5 - (-12) = 17 
b ) | "4 _g| = 27 - 8 = 19 
C) | J "5| = 35 - 18 = 17 
9.02 a) | g'g 3 | = 8,4 - 16,8 = -8,4 

b) | 4^2 ~ l ',0 | = 13 ' 35 + 13 ' 44 = 26 ' 79 

c ) | ~ o ' f 17 | = °' 1224 + °' 1162 = °' 2386 


I 3 7 -2I 

9.03 a) 8 4 5 = 3-4-6 + 7-5-(-2) + (-2)-8-(-3) 

| -2 -3 61 

- ( (-2)■4•(-2) + 3•5•(—3) + - 7-8-6 ) = -257 


-4 3 2 

b) -5 2 7 = 24 + 168 - 60 - 32 + 168 - 45 = 223 

8 6-3 


5-18 

C) 6-3 4 = -75 -16+96+96-40+30 

4 2 5 = 222 - 131 = 91 



a i2 3 11 
a 22 a 21 


Zu 3 


a b 
a b 


a b - b a 
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9. Determinantenrechnung 


a ll k a l 2 
a 21 k a 22 


- a l 2 a 2l> = k 


ka ll ka l 2 
a 21 22 


"11 "22 12 21 


; k a.. a__ - k a. _ a, 


12 21 


= k < a ll a 22 " a l2 a 21> * * 

| = k < a il a 22 - a i2 a 21> 

I k 

(k a il> a 22 - < k a i2 } a 21 " 


a il a i 2 
a 21 a 22 


“11 “12 
a 21 a 22 


11 12 
a 21 a 22 


"11 "12 
k a il k a i 2 


: k a ll a l 2 - k a l 2 a il = 0 

für jede Belegung von a ni und a r 


9.05 a) 


2 

6 


5 

7 



10 - 7 - 4-(4-6) + 3-(14 - 30) 

3 + 4-2 - 3-16 = 3 + 8 - 48 = -37 


1 1 4 3 1 

2 5 1 = 1-5-2 + 4-1-6 + 3-2-7 

I 6 7 2 I 

-(3-5-6 + 4-2-2 + 1-1-7) 


= 10 + 24 + 42 - (90 +16+7) 
= 76 - 113 = -37 


14 3 

2 5 1 = 

6 7 2 


0 -3 -5 

0 -17 -16 




9. Determinantenrechnung 


37 


= 414 - 152 = 262 


= 5 (12 - 18) - 6 (-9 + 16) + (-27 + 32) 


= -30 - 42 + 5 = -67 


= -140 - 24 + 72 - (90 - 64 - 42) = -92 + 16 


= 380 - 64 = 316 


- (-100 - 336 + 32) 


3 0,24 0,8I I 3 24 8I I 3 0 23I 

C) -2 6,7 1 = 0,001- -2 670 10 = 0,001-2 686 0 

| 5 0,2 4,3| | 5 20 43| | 5 -20 68| 

= 0,001-[686 (204 - 115) +20 (0 + 46)] 

= 0,001-(686-89 + 920) = 61,054 + 0,920 = 61,974 
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9. Determinantenrechnung 


I 3 0,24 0,8I 

-2 6,7 1 = 12,9-6,7 + 1,2 - 0,32 - 26,8 + 0,48-4,3 

| 5 0,2 4 , 3 [ 

= 86,43 + 1,2 - 27,72 + 2,064 = 61,974 


9.07 a) 



y 1 (x 3 - x 2 ) + y 2 (-x 3 + x x ) 
X 1 (y 2 " y 3 ) + x 2 (y 3 " y l } 


y 3 (x 2 - x 3 ) 
x 3 (y l" y 2 } 



c) 


5 9 7 
2-9 4 

8 7 -24 


7 0 11 7 0 11 

2 -9 4 = 2 9 4 

8 7 -24 0 43 40 


2520 + 946 - 1204 = 3466 - 1204 


2262 


15 9 71 

2 -9 4 = 45-24 + 36-8 + 98 + 63-8 - 140 + 18-24 

|8 7 —24| 

= 63-24 + 99-8 - 42 = 1512 + 792 - 42 = 2262 




9. Determinantenrechnung 
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9. Determinantenrechnung 



= -3 (32 + 36 - 15 - 48 + 12 - 30 ) =39 

I 4 6 l| . - 

= -3 1 2 0 = -3 f = -3 (-9 - 4) = 39 

2-9 0 I 2 9 I 


C) 5 0 0 0 

9 q 9 q ist eine Diagonaldeterminante. 

0 0 0 10 

Der Wert einer Diagonaldeterminante ist gleich dem Produkt 
der Elemente der Hauptdiagonale: 

D = 200 



10. Matrizenrechnung 


10.01 a) 30 b) -108 c) 120 

10.02 Entwickelt man die Dreiecksdeterminante nach der ersten 

Spalte, erhält man a^- a.^, da alle weiteren Produkte Null 
werden, kann nun wieder nach ihrer ersten Spalte ent¬ 

wickelt werden. Dies führt auf a ^' a 2 2 2 a 22 und schließlich 
aUf a ll a 22 a 33 '’’ a nn 

(Hier wurde die Adjunkte zum Element a 22 der adjungierten 
Matrix mit 2 a 22 beze; *- clinet *) 


10.04 a) -2 

10.05 Für: = 1 a 



0 0 

0 0 ) 00 

0 0 J 0 0 

0 0 

b) 43 



0 0 0 0 ) 
0 0 0 0 J 


C) 


-217 


a 21 


6 



0 

6 


C T = (3 4 -5 2) 


10.08 r = 1, wenn a 13 = a 23 a (a lt * 0 v a 12 * 0 v a 13 * 0) 

r = 2, wenn a 13 * a 23 a (a 1;L * 0 v a 12 * 0) 

10.09 Es liegt folgendes lineares Gleichungssystem vor: 

a ll X 1 + a 12 x 2 + . + a ln x n = b l 

a 21 X 1 + a 22 x 2 + . + a 2n x n = b 2 


a x. + a x,+.+ a m „ x„ = b 

ml 1 m2 2 mn n m 

Die Koeffizientenmatrix lautet: 

A = (a. k ) i-1, 2, m 

k = 1, 2, ..., n 
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10. Matrizenrechnung 


Die erweiterte Matrix A g erhält man durch Hinzufügen 

des Spaltenvektors ( t> i ) in die Koeffizientenmatrix. 

Rang A = Rang A g = n: Es gibt ein Lösungs-n-Tupel. 

Rang A = Rang A g < n: Es gibt ® viele Lösungen. 

Rang A < Rang A g : Die Lösungsmenge ist leer. 

(Widerspruch im Gleichungssystem) 

Rang A > Rang A g : Dieser Fall kann nicht eintreten. 

10.10 Rang A = 2 Rang B = 3 Rang C = 3 

10.11 Hier gibt es viele Möglichkeiten von Beispielen, besondere 

herauszustellen erscheint hier nicht notwendig. 

1. Folgt aus der Definition des Ranges einer Matrix. 

2. Eine Determinante gibt es nur für eine quadratische 
Teilmatrix. 

3. Dies kann nur eine Änderung des Vorzeichens einer für den 
Rang maßgebenden Determinante bewirken. 

4. Dies bedeutet für eine, für den Rang maßgebliche Deter¬ 
minante ebenfalls nur eine solche Vertauschung von Spal¬ 
ten und Zeilen. 

5. Durch diese Änderung kann der Wert der Entscheidungs- 
determinante nicht 0 werden. 

6. Siehe 5. 

7. Folgt aus den Regeln für das Berechnen des Wertes einer 
Determinante. 

10.12 a) [ 6 2 4 1 b) ( 4 -8 0 ) e) ( 8 -3 2 | 

3 4 4 -1 4 -2 1 7 1 

[l3 19 5 J l-l -9 -19 J [ 6 5 -4 J 

c) (5 3 17 5 8 ) d) f-3 -1 7-1 -6 ] 

U -1 -4 1 4 J l 2 -1 8 5 -4 J 

10.13 a) [ 13 53 25 ) b) ( 22 27 -7 ) c) f 79 29 ) 

16 19 26 28 9 -17 [-71 -15 J 

[-33 -68 -57 J [l21 95 -56 J 

10.14 a) { 15 -37 -2 ) b) f 22 10 6 14 40 ) 

-6 16 -11 [ 0 2 -40 -18 24 J 

[-11 -50 -88 J 

c) f 2 -9 -19 ) 

-14 4 -41 

[ 0 -14 -43 J 




10. Matrizenrechnung 
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f 5 -25 0' 

c) -10 30 0 

- 5 9-4 


0 H 0 

10.18 a) 28 9 -35 

" 21 -7 -7 


12 -6 -14 
6 3-5 

0 0 8 


. 0 0 4 

c) 8-4 -13 

0 8 10 


Vergleich mit 


18 24 -34 

3 -12 23 


[ -6 24 -38 

= 4I -18 24 _34 

40 3 -12 23 


=> (B a -A - a ) • (A-B) = E | . (A-B) 


(B -A ) • (A-B) • (A-B) " 


E- (A-B) =* B -A = (A-B) 
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10. Matrizenrechnung 




11. Integrationsmethoden 
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11. Integrationsmethoden 




11. Integrationsmethoden 


b) I = (x* + 2 x) ✓ 1 + x“ dx 


(sinh t + 2 sinh t) cosh t dt 


= sinh t cosh t dt + 2 cosh t sinh t dt 


I. = sinh t cosh t dt ( siehe: 11.03 a) mit a=l ) 


* + 1 (2 x J + x) - ln|x + V x*+ 1| 


cosh t sinh t dt 


2 I udu = 2 ^- + C_ 


ln|x + v 1 + x 


a sinh t dt + a 


t _ va + x + x 


du 

J u 2 - 1 


a-cosh t 


a ln 1 + e . + C 



11. Integrationsmethoden 


[ dx = f 2 dt 
J cos x J 1 - t 2 


✓l + cos x + VI - cos x 
✓1 + cos x - VI - cos x ' 


1 + sin x 
cos x 


sin x cos x J sin 2 x 


ln|tan x| + C 


f sin 2 x + i 
J sinx ca 


j(tan x + cot x) dx = ln|sin x| - ln|cosx| + C 
= ln|tanx| + C 


Hier führt tan x = t schnell zum Ziel! 


= 1 ^ 


sin x 





11. Irttegrationsmethoden 


oder (unter Verwendung von 11.04 a)): 

f sin 2 x_+_ cos^_x dx = f _d^_ + f j 
J sin 2 x cos x J cos x J 


|1 + sin x| 1 


11.05 a) INT(l/(SIN(x) . COS(x)) A 2, x) 


b) INT(1/(2 + 5 COS(x)), x) 


r-LN r. (^ :. 9? s( x) ,_- 

L(V r 7COS(x) + 


✓3-SIN(x) + Vl~\ 


(X) + ✓3-SIN(x) + V7-I 


C) INT(1/(2 SIN(X) + 3 COS(x)), x) 

VI3 • LN [K13 - 2) • COS (x) + 3 • SIN (x) + /13 - 21 
L(v' 13 + 2) - COS(x) - 3-SIN(x) + /13 + 2J 


1.06 a) I n = | x n e x dx = | x n d(e x ) = x n e x - j e x n x"' 1 


b) I„ = | cos 11 x dx = | cos 11-1 x < 


= sin x-cos x + (n - 1) 


x (1 - cos x) dx 


I = sin x-cos x + (n - 1) I _ - (n - 1) I 


I 1 = sin x + C I Q = x + C 

I n = | x n sin x dx = - | x n d(cos 

= - x" =o. X + „ | 


n-1 


cos x dx 
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11. Integrationsmethoden 


I = - X COS X 


• n | x"- 1 , 


I = - x cos x + n x sin x ■ 
I, = - x cos x + sin x + C 


<1) I n - | x" • 


= x sin x + cos x + C 


dx = | ln|x - 2| 


c) 

11.08 a) 

b) 


a) f 3 x - 5 

' J (x - 2) (x + 4) 

b) j (X + 5) (X - 2) dx = — ln l x + 5 
f 4 X - 5 

J X 2 + 2 X - 

f 10 X + l 

J x 2 - 7 x + 

I 


■ dx = ln|x + 5| 


12 


3 x - 2 x + 1 


■ dx = 45 ln|x - 4| 


x (x - 5) (x + 7) 

ln | x — 5 | + yk ln | x + 7 | 


10 


dx 


f 4 x + 5 

J (x - 1) (x + 1) (x - 6) 

= _g. lnl x ~ ll + 29 ln|x - 6| + 
10 35 

[ 2 X ~ \ dx = —-— + 2 ln|x + 

J (x + l) z x + 1 


X) 

n (n - !) I n _ 2 
: Q = - cos x + c 

(n - 1) I n _ 2 

I Q = sin x + C 
+-^ln|x+4| +C 
i| + |ln|x - 2| + C 
+ | ln|x - 3| + c 

- 35 ln|x - 3| + C 

' 3 § ln|x| + C 

ln|x + 1| 

14 

1| + c 


+ c 
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= 5 ln|x - 1| 


x - 1 (x - 1)* 


C, [ 3_jd..--y + 1 dx = ln M + 26 ln l * - 3 I-25_ + c 

J X (X - 3) 2 9 9 3 (X - 3) 


11.10 a) [ - x!jL_ 1- dx 

J X (X - l) (X + 2) 

_ _1_ 7 ln|x + 2| 2 ln|x - l| _ ln|x| 

2 x 12 + 3 4 u 

*» 1 = [ dx 

J x J + X 4 

X 2 + 1_A,B .C D 

3 4 + 2 + 3 + 

X + X X X X X + 1 


x 2 + 1 = A (x 3 + x 2 ) + B (x 2 + x) + C 
Koeffizientenvergleich: A + D = 0 
B + C = 0 C = 1 B = - 1 


I 



X + 1 

2 x - 1 


(x + 1) + D x 3 
A + B = 1 
A = 2 D = - 


2 


c) I 


[ 3 x 2 - X + 2 

J X 2 (X + l ) 2 


dx 


3 x 2 - x + 2 _ A + B_ + C + D 
x 2 (x + l) 2 x x 2 x + 1 (x + l) 2 

3 x 2 - x + 2 = 

A (x + l) 2 x + B (x + l) 2 + C (x + 1) x 2 + D x 2 
x = 0 => B = 2 x = -l^D=6 

x=l=>2A + C = -5 x=2=»3A + 2C = -5 

C = 5 A = - 5 
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11 . Integrationsmethoden 



x + 1 (x + l)“ 
_ 8 x + 2 * „ 


dx 


Lösung mit Hilfe von MATHEMATICA 

a) Integrate[(x A 3 + 1) / (x A 2 (x - 1) (x + 2)), x] 

1 2 Log[l - x] Log[x] 7 Log[2 + x] 

-+---+- 

2X3 4 12 

b) Integrate[(x A 2 + 1) / (x A 3 + x A 4), x] 

-1 1 

- 5 + - + 2 Log[x] - 2 Log[l + x] 

2 X^ X 

c) Integrate[(3x A 2 - x + 2) / (x A 2 (x + 1) A 2), x] 

-2 6 

- 5 Log[x] + 5 Log[1 + x] 

X 1 + X 


Partialbruchzerlegung mit Hilfe von MATHEMATICA 

a) Apart[ (x A 3 + 1) / (x A 2 (x - 1) (x + 2)) ] 

2 117 

- -+- 

3 (-1 + x) 2 x 4 x 12 (2 + x) 

b) Apart[ (x A 2 + 1) / (x A 3 + x A 4) ] 

112 2 

x 3 x 2 x 1 + x 

c) Apart [ ( 3x A 2 - x + 2) / (x A 2 (x + 1) A 2) ] 

2 5 6 5 

2 X 2 1 + X 

X (1 + x) 


11.11 a) I 


I 


X + 2 

x 2 - 4 x + 7 


dx 


I 


1 r 2*- 4 dx + 4 [ *4— 

2 J x^-4x + 7 J (x - 2 p + 3 


= - ln|x 2 - 4 x + 7| + I. 
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Z 4 [ t_ _ 4_ r dt_ 
1 J 3 t 2 + 3 VI J t 2 + 1 


= — arctan x ~ 2 

VI VI 


x = 2 + VI t 


dx = VI dt 


I =-ln|x 2 -4x+7|+— arctan - 


b) I 


, f 3 X - 5 

J x 2 + 6 x + 


(2 X + 6) • 1,5 _ 14 


c) I 


x + 6 x + 12 x + 6 x + 12 (x + 3) + 3 

= 3 . 2 x + 6 _14_ 

2 X 2 + 6 X + 12 3 (1 + t 2 ) 

! = 2 [ ** + « - dx _ 14VI f _dt 

2 J xS 6 x + 12 3 J 1 + 

= - ln|x 2 + 6 x + 

2 

, [ 1 X - 5 ( 

J x 2 - 2 X + 4 


X + 3 = VI t 


ln Ix + 6 x + 12 I-arctan - 

VI VI 


= — lnlx - 2x + 4| + — arctan —-— + C 

2 VI VI 

1.12 a) _ 8 x ~ 1 dx = 4 ln| x 2 + x + 2 | + — arctan 
J x z + x + 2 V7 

b ’ 1 = 1 * -1 


2x + 1 
V7 


x + 1 


(X - 1) (X + X + 1) 
A_ , B x + C 


(X - 1) (x“ + X + 1) X - 1 X + X + 1 

x 2 + 1 = A (x 2 + x + 1) + B (x 2 - x) + C (x - 1) 
x = 0; 1; -1 liefert: A = 2 B = 2 C = 
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11. Integrationsmethoden 


I = 


x - 1 


3 (x - 1) 3 (x+x + l ' 

I = — ln|x - 11 + — f j X ~ 1 dx 

3 3 J X 2 + X + 1 

2 ln|x - l| + ln|x 2 + x + l| 


[ x 2 + 1 

J x 3 + 1 


2 ln|x + 1| + ln|x 2 - x + 1| + 


a) [ l* 2 2 " ' " dx = f (- V ~ 1 - 

J x 4 + x 2 + 1 J ^2 (x 2 - x + 


1 arctan 2x v/3 ~ + < 


1) 2 (x + x + 1) J 


x - x + 1 

' ln ~2 - 

1 x+x + l 1 

— (arctan 2x + 1 + arctan 2x ~ 1 
21/3" v V? 


Lösung mit DERIVE 

INT( (2 x A 2 - 1) / (x A 4 + x A 2 + 1), X ) 

✓5 atan[^-^^I ] V3 atan[ ^ x + i) 


3 ln r x l ~ x +1 1 

L x + x + 1 -I 


Partialbruchzerlegung mit Hilfe von DERIVE 
EXPAND(2 x A 2 - 1) / (x A 4 + x A 2 + 1) 

_3 _x_1_ 3 x 

2 (x 2 - x + 1) 2 (X 2 - X + 1) 2 (X 2 + X + 1) 


2 (x 2 + x + 1) 
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b) f 2 V + 3 

J x z (x z - X + 


1) 


= - ln ■ 


• arctan - 


x + 1 x VS VS 

Lösung mit DERIVE 

INT( (4 x + 3) / (x A 2 (x A 2 - x + 1)), x) 


VS ATAN 


p? (2 x - 1) 1 

L_ 3 -I _ 7 LN(x z - x + 1) + ? . 


Partialbruchzerlegung mit Hilfe von DERIVE 
EXPAND( (4 x + 3) / (x A 2 (x A 2 - x + 1)) ) 

_Z_J<_ + _4_ + _3_ + 1 

2 2 2 
x-x + 1 x-x + 1 X X 


= [ 2 « - 5 * 7 d» 

J X (X + X + 1) 


2 x + x -5 
x (x 2 + x + 1) 


x + x + 1 


2 x + x - 5 =A (x + x + 1) + B x + CX 
x = 0 , 1,-1 liefert: A=-5 B = 7 C = 6 


- f (- + l xtt ] *> 

J L X x z + x + V 
: = -5 ln|x| + [ 7 x + 6 dx 

J x + x + 1 

= -5 ln|x| + 3,5 ln(x 2 + x + 1) 


— arctan - 

VS VS 


Lösung mit DERIVE 

INT( (2 x A 2 + x - 5) / (x (x A 2 + x + 1)), x ) 

5 ATAN ^' ( 2 x + !) ] 2 

L 3 J . 7 LN (x 2 + x + 1) 

3 2 5 

Partialbruchzerlegung mit Hilfe von DERIVE 
EXPAND( (2 X A 2 + x - 5) / (x (x A 2 + X + 1) ) 

7 X 6 5 

2 + 2 “ 
x+x + 5 x+x+5 x 



12. Unendliche Reihen 


12.01 a) Die Reihe ist als Majorante zur harmonischen Reihe 
divergent . 

b) wie a). 

c) Zunächst: Die Reihe + yy- + ...ist konvergent. 

1 = 1 1 
n (n + 1) n"n + l 


= 1 - n-T-T 

s = lim s = 1 
n -> oo 

Für unsere Reihe - ohne das erste Glied - ist die obige 
konvergente Reihe Majorante. 

Daher konvergiert unsere Reihe. 

12.02 a) Diese Reihe konvergiert , weil die folgende geometrische 

Reihe Majorante ist: 1 + - + + ... (q = - < 1) 

5 5^ 5 

v n € IN*: —— ^ 

5 n n 5 n 

b) Diese Reihe divergiert , weil sie Majorante zur har¬ 
monischen Reihe - ohne das erste Glied - ist : 


c) a fi = n n — Von der gegebenen Reihe ist jedes Glied 

größer als das der harmonischen Reihe 1 + j + y + ..., 
die daher eine divergente Minorante ist. 

Die gegebene Reihe ist daher divergent. 

12.03 a) Diese Reihe konvergiert , weil : 

a) lim --- = 0 

n ■> « 2 n - 1 

• i i i i i i 

ß) V n e IN 


2 n 
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b) 


Diese Reihe konvergiert , weil : 

a) lim --- ^ = 0 

n m (2 n - lj 


ß) wie oben 


c) Diese Reihe konvergiert , weil : 

a n 

a) < 1 — - > für a < 1 eine Nullfolge ist und 
• | a n | | a n+1 | 

ß) V n e W ; |- > | n + 2 | gilt, weil das 

Folgeglied sicher keinen größeren Zähler, aber 
einen größeren Nenner hat. 

d) Das Kriterium von Leibniz liefert keine Entscheidung , 
weil : 

V n e IN*: 2 — + ^ (Man kann die Divergenz be¬ 

weisen) 

e) a) Der Betrag des Folgegliedes wird hier durch Multipli¬ 

kation mit dem Faktor ( 2 ' n - 1)^2 n - 2) gebildet. 
Dieser ist für n fc 2 sicher kleiner als 0,7; daher 
ist lim u = 0. 
n -> co n 

«I l»nl - l u n.ll ‘ * (2 n ? 1)-2 n ) » 0 

für n s 2 

daher ist die Reihe konvergent . 


12.04 a) (1) Die konvergente Vergleichsreihe (siehe: 12.01 c)) ist 
Majorante für die Beträge der gegebenen Reihe für 
alle Werte von x, da die Kosinusfunktion (im Reellen) 
nirgends größer als 1 wird; die Reihe konvergiert 
absolut für alle x e IR . 

|sin n(x + n/4)| 

(2 ’ '“n 1 -^^ 

n n 
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12. Unendliche Reihen 


Die Reihe 1 + ^ + ... konvergiert (Aufgabe 

7T 3^ 

12.01c) und ist Majorante der vorliegenden Reihe. 
Daher konvergiert die vorliegende Reihe absolut. 


12.05 .) lim = 1 


*» “■ (ÄV 2 )! n - 


Keine Entscheidung mit Quotienten¬ 
kriterium (Divergenz ist zu beweisen). 


c) lim 


n + 1 


(n + 1)! ‘ nT 
Die Reihe ist konvergent . 

•(2 n - 1 ) 
. • (2 n) 


= lim - = 0 

n->oo 


a > iäQKSH TT 


= lim 
n-»® 


2*4*6* ... 
[ 2 n - 1 


•< 2 » - 3 ) . 1 )] 

. •(2n - 2) n JJ 


1 = 1 


2 n n + 1 I rt 2 n + 2 

' n->w 

Keine Entscheidung m it dem Quotientenkriterium. 


12.06 a) R„ +1 - (*^i * 2 " +1 0 = « a 1 

Beachten Sie: Das n-te Glied hat die Form 


. . „ cosh (i> x) 2n 

b > R n+1 = (2 n)1 x 


0 s s 1 


(2 n)! 

Beachten Sie: Das n-te Glied hat die Form 
,.n-k k 


c) R k+1 ' 


(1 + x) 1 


0 s s l 


x 2n_1 
(2 n - 1) ! 


x 2n-2 
(2 n - 2) ! 


12.07 a) f c jj |^2 + 8 tan a + 6 tan aj ( = 1 ) 

a 2 |l + 4 a 2 + 3 a 4 j = 0,003 mit tan a = a 

a < 3,14° ( und nach Newton verbessert a < 3,12° ) 
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b > f < h 

---=- < 0,001 

- (-3 

a < 22,12° 

a) a < 14,0° 

b) a < 9,4° 

c) a < 37,8' 

a) 8,26° 

b) 3,83° 

c) 1,78° 

a) 32,6° 

b) 20,5° 

c) 13,0° 

1,67 • IO -8 = 0,0000000167 


sinh(jx) = j 

„ + -s3 x 3 .5 x 5 

x+D 31+] 5l + -** 



f x 3 x 5 1 


= j 

[ X " 31 + 5T ‘ * * * J 


= j 

sin x 


cosh(jx) = 1 

, •2 x .4 x 

+ 1 JT + 1 4 T + 


= 1 

x 2 x 4 

' JT + 41 " 



cos x 



13. Fourier-Reihen 


2tr 

0 J 


für 

für 


r 21 

■sj- 

TT" 


f(x) cos nx dx 


r 2n 

b n - I | f< * : 

_ a 
n 

b „- { 


■*f 

n 


f(x + 2n) 


cos nx dx : 


sin nx dx 


1 2n 


‘ is U-» n - l] 


-2a 
n n 


f(x) = f ■ 


für ungerade n 
0 für gerade n 

2a f sin x sin 3x sin 5x 


13.02 1. Möglichkeit: Vergleich mit Aufgabe 13.01: 

a Q bleibt, der Graph ist an y = a Q gespiegelt, 
daher Vorzeichen der Sinusglieder umkehren und 
wegen der Halbierung der Periodenlänge x durch : 
ersetzen: 

_ 2_a 

2 


f(x) = f + 


sin 2x 
1 


sin 6x 
3 


2. Möglichkeit: Wenn die primitive Periode n ist 
ist 2 n auch eine Periode: 


f(x) = | 

f(x ± 2TT) 


a für 0<x<- 
0 für j < x < i 

= f (x) 


und für TT < x < —^ 
und für < x < 2 tt 


f (x) 


., dann 


a = a 
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a-sin nx dx + 


a- sin nx dx 


cos nx + cos nx 

n lir/2 


' 1 3TT/2 

+ cos nir - cos - 


f(x) 


(-1) - 2-cos n7i • cos — j- J 

■ JTÜ [ 1 + (- 1 )" * 2 (- 1 )" cos V 1 ) 

für n = 4 1 + 2 (i e N Q ) 

für n = 2 1 + 1 und für n = 4 1 
2a f sin 2x , sin 6x . sin lOx . ) 


( ±_a 

-1 r 


13.03 Dies ist Sonderfall 2. von Beispiel 4. 

13.04 Gerade Funktion: f(x) = x für 0 < x < n 

f (x) =2n - x für 7r < x < 2n 
f(x) = f(X ± 271) 
f 2 77 pTI 2 . TT 

a o " S J f < x > dx = iJ xdx = |r^| =Tr 


s-if 


x cos nx dx = 


.277 

I ( c ' 


■ x cos nx dx 


x sin nx , cos nx 

-Z- _2 


x sin nx cos nx 


(- 1 )“ - 1 + (- 1 )^ 




n n 

cos x cos 3x cos 5x 


für gerade n 
für ungerade n 



62 


13. Fourier-Reihen 


13.05 Ungerade Funktion: a Q = a n = 0 n e IN 

f(x) = a - für 0 < x < 2 tt und f(x ± 2n) = f(x) 


.271 

b n = I J (nr^jsinnxdx 


n x cos nx - sin nx 


13.06 Ungerade Funktion: a Q = a n = 0 (n e IN) 

4 a x 

I 

f (x) 


(TT - X) 


4 a x 


für 0 < x < tt 

für - 7T < x < 0 
n“ 

f(x ± 2 tt) = f(x) 

i r" 2 

b n = — f(x) sin nx dx = — f(x) sin nx dx 
—71 0 

= | (tt x - x 2 ) sin 


sin nx dx = 


(sin nx x cos nx 




cos(nx - n 2 x 2 ) + 2 n x sin nxj 
für gerade n 


32 a 
3 3 


für ungerade n 


_ 32 a ( sin x + sin 3x + sin 5x + 
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13.07 Gerade Funktion: f(x) = j für -n s x s u 
f(x) = f(x ± k- 2n) ] 

»0 - s { f <*> a * - I } i a * - Hr 
0 0 

r 71 2 

a n - I j I cos nx dx - i J 2-f cos mrj - ( 
0 

. _ TT 2 ( COS X COS 2x COS 3X 


Ungerade Funktion: f(-x) = - f(x) 

f(x ± k- 2n) = f(x) 
a. = a_ = 0 (n e IN) 


b n - I [ 


f(x) sin nx dx 


p7i-a tt- a 

sin nx dx + b sin nx dx + ^ (x-tt) sin nx dx 
^ a •’n 


J x sin n x dx = 


x cos nx , sin nx 


j; 


x sin n x dx = 


a cos na sin na 


sin nx dx = cos na - cos n(n-a) 


J (x - tt) sin nx dx-ÜL.-.»> cos nx + ^in_nx 


r -a 

J (x ■ 71 ) 


sin nx dx = • 


i cos n(n-a) + (~l) n sin i 


/ .vn+1 . 

sin na + —'—= sin na 
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13. Fourier-Reihen 


Es ist somit: 


2 b 
a Ti n 2 


0, falls n gerade 
4 b sin na für un- 
a tt 2 gerade 

n 


- i-t ( ^ si " + •••) 


Sonderfälle: 


4 b f sin a 


Formal erhalten wir: 
. ,sin 3a , 


a 0 =» lim - 

a-»0 


Die Fourier-Reihe der Dreieckskurve lautet: 
... _ 8b ( sin x , sin 3x sin 5x ) 


13.09 In der Hilfsdatei INT_APPS.MTH ist die Funktion 
FOURIER(y, t, tl, t2, n) definiert. 

F0URIER( ABS(SIN(x)), x, -tt, tt, 3) 

COS(3 x) _ 5 COS(2 x) 2 

4 TT 3 TT TT 

Wir zeichnen nun die Kurven mit der VECTOR-Funktion. 
VECTOR( FOURIER( ABS(SIN(x)), x, -n, TT, n ), n, 3 ) 

In allen drei Fenstern sehen wir zum Vergleich y = |sin x| 
Zusätzlich im Fenster 
1 den Graphen von ^ 

3 den Graphen von ^ 

4 den Graphen von § - 5 ^ x ) + C °|% *) 




Nun lösen wir das Beispiel mit MATHEMATICA. 
Needs["Calculus'FourierTransform'"] 
f[x_]:= Abs[Sin[x]]; 

gl=Plot[f[x]/.x->Mod[u,2Pi],{u,-2Pi,3Pi}] 



n2=Normal[FourierTrigSeries[f[x],{x,0,Pi},9 ]]//N 
g2=Plot[n2,{x,-2Pi,3Pi}] 

0.63662 - 0.424413 Cos[2. x] - 0.0848826 Cos[4. 
0.0363783 Cos[6. x] - 0.0202102 Cos[8. x] 

- 0.012861 Cos[10. X] - 0.00890377 Cos[12. x] 

- 0.00652943 Cos[14. x] - 0.0049931 Cos[16. x] 

- 0.00394192 Cos[18. x] + ... 







14. Extremwerte bei Funktionen 
in zwei Variablen 


14.01 Wir lösen die Aufgabe zunächst mit Hilfe von DERIVE. 
Z(x»y):=3 x A 2 - 2 y A 2 +6x+4y+12 
DIF(3 x A 2 - 2 y A 2 + 6 x + 4 y + 12, x) 6 X 

DIF(3 x A 2 - 2 y A 2 + 6 x + 4 y + 12, y) 4 - 

SOLVE([6 x + 6 = 0, 4 - 4 y = 0], [x, y]) [ [X = -1, 

Z(-l, 1) 11 

Nun prüfen wir mit der hinreichenden Bedingung 

.. , _ ä2 < x 0'V s2(I< 0' y 0 ) f ä2 < x 0' y 0> l 2 

4(« 0 .y 0 ) ^ (“ x - ay J 

ob ein Extremwert vorliegt. 



DIF(6 x + 6, x) DIF(4 -4 y, y) - DIF(4 -4 y, x) A 2 
Es liegt an der Stelle (—1|1|11) sicher kein Extremum 
f-llllll) ist ein SattelDunkt. 


Leichter geht es mit MATHEMATICA. 


Hier steht uns die Funktion FindMinimum zur Verfügung. 


z[x_, y_] := 3 x A 2 -2 y A 2 + 6 x + 4 y + 12; 
FindMinimum[z[x, y],{x, -l},{y, 1}] 

FindMinimum::fmgz: 

Warning: FindMinimum encountered a vanishing gradient 
The result returned may not be a minimum; it may be a 
maximum or a saddle point. 





14. Extremwerte bei Funktionen in zwei Variablen 





{11., {x -> -1., y -> l.}} 

Plot3D[z[x, y],{x, -2, 1}, {y, -1, 3}] 


14.02 Minimum (-2|1.333|-3.333) 
z[x_, y_] := 

X A 2 + 3 y A 2 + 4x-8y+6; 
FindMinimum[z[x, y], 

{x, -l}, {y, l}] 

{-3.33333, 

{X -> -2., y -> 1.33333}} -4 

Plot3D[z[x, y], {x, -4, 0}, 

{y, -i, 3}] 


14.04 Minimum (-5|0.0833|-30.0833) 
z[x_,y_] := 

x A 2 + 12 y A 2 - 2 y + 10 x - 5; 

FindMinimum[z[x, y], 

{x, -4}, {y, 1}] 

{-30.0833, 

{X -> -5., y -> 0.0833333}} 


Plot3D[z[x, y],{x, -20, 20}, 
{y, -20, 20}] 


14.03 Sattelpunkt (0|0|10) 
z[x_,y_] := 

x A 2 + y A 2 -3 x y + 10; 
FindMinimum[z[x, y], 

{X, 0}, {y, 0}] 

{10., {x -> 0., y -> 0.}} 


Plot3D[z[x, y], {x, -15, 5} 

{y, - 2 , 2 }] 
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14. Extremwerte bei Funktionen in zwei Variablen 


14.05 Wir lösen die Aufgabe zunächst mit Hilfe von DERIVE. 

3 ~ x A 2 / 16 - y A 2 / 9 

DIF(3 - X A 2 / 16 - y A 2 / 9, x) -0.125 X 

DIF(3 - x A 2 / 16 - y A 2 / 9, y) -2/9 y 

Leicht ist zu erkennen: (o|o|3) ist der gesuchte Punkt. 


Aus der Zeichnung ersehen wir: 
(0|0|3) ist ein Maximum 
Wir überprüfen: 

q2 

2-| = - 0,125 < o Maximum 
dx 



Nun mit MATHEMATICA: 

*[*_» y_] := 

3 - X A 2 / 16 - y A 2 / 9; 
FindMinimum[-z[x, y], {x, 0.5}, 
{y, o.5>] 

{-3., {X -> -1.25455 10 _1 ? 

y -> -2.67134 10 -13 }} 

Plot3D[z[x, y], {x, -2, 2}, 

iY, - 2 , 2 }] 



14.06 Sattelpunkt(0|0|0) 

z[*_, y_] := 

(2 y A 2 - x A 2) / 4j 



FindMinimum[z[x, y], 

{*, 0}, {y, 0}] 

FindMinimum::fmgz: 

Warning: FindMinimum .... 
it may be a maximum or a saddle 
point. 

{0., {x -> 0., y -> 0.}} 

Plot3D[z[x, y], {x, -2, 2}, {y, -2, 2}] 



14. Extremwerte bei Funktionen in zwei Variablen 


14.07 Maximum(4|-3|-2) 

z[x_, y_] j = 2 - (X - 4) A 2 / 16 - (y + 3) A 2 / 9; 
FindMinimum[-z[x, y], {x, 0>, {y, 0>] 

{-2., {x -> 4., y -> -3.}} 

Plot3D[z[x, y], {x, 2, 6}, {y, -5, -1}] 



14.08 E(± n n/6 | ± n tt/4 | ± 2) 

z[x_,y_] := 2 Sin[3 x] Cos[4 y]j 
FindMinimum[z[x, y], {x, 0>, {y, 0}] 

{-2., {X -> -0.523599, y -> 0.}> 

Flot3D[z[x, y], {x, -Pi/2, Pi/2}, {y, -Pi/2, Pi/2>] 


Nun lösen wir dieses Beispiel mit DERIVE. 




16. Differentialgleichungen erster Ordnung 


16.01 a) y = K e 

c) y = ± /2 e x + C 


16.02 a) y = / 3 e x (l - x) + K 
c) y = 2 (3 x 2 + K) _1 

16.03 a) y = C 1 (x - 2) - 1 

2 y/yr 

c) y = 0,541 e* vx 


b) y = ± / K - 


b) y = 6 e 


c) y = 


= - 2 e 1 ' 2 - 2 / x 


16.05 1. DSOLVEl_GEN(-y, 1, X, y, c) 

2. DS0LVE1_GEN(- x y, 1, x, y, c) 
DSOLVE1(- x y, 1, x, y, 0, 2) 

3. DSOLVE_GEN(x, y, x, y, c) 

4. Im Buch (Seite 118) gelöst. 


y = k e 

x 2 

LN(y) - | = c 

y = 2 exp(f 2 ) 

x 2 y 2 
$ + i = c 


5. DSOLVE_GEN(-l - y, x, x, y, c) y = k x - 1 

6. DSOLVE_GEN(- k y/x, 1, x, y, c) y = C x k 

7. DS0LVE1(- 1 - y 2 , 1, x, y, p/4, r3) y = TAN(x + ^ ) 

8. DSOLVE(1 n, 1, t, n, 0, nO) 

LN(n) - LN(nO) + A n = 0 

SOLVE(LN(n) - LN(nO) + 1 n = 0; n) 

( n = nO • e” X 

9. DS0LVE1(q/c, r, t, q, 0, p) 

- r LN (p) + r LN(q) + - = 0 q = p e _t/(c r) 
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16.06* l. DSolve[y’[x] == y[x], y[x], x] 

{{y[x] -> E A x * C[l]}} 

2. DSolve[y*[x] == x y[x], y[x], x] 

{{y[x] -> E-(x-2/2)*C[l]}} 

"Anfangsbedingung y(0) = 2:" 

DSolve[{y*[x] == x y[x], y[0] == 2}, y[x], x] 

{{y[x] -> 2 * E-(x-2/2)}} 

3. DSolve[y[x] y’[x] + x == 0, y[x], x] 

{{y[x] -> -Sqrt[-x 2 + C[l]]>, 

{y[x] -> Sqrt[-x 2 + C[l]]}> 

4. DSolve[y[x] y’[x] + 2 == 0 , y[x], x] 

{{y[x] -> -Sqrt[-4 x + C[l]]}, 

iy[x] -> Sqrt[-4 x + C[l]]}} 

5. DSolvefx y’[x] == 1 + y[x], y[x], x] 

{{y[x] -> -1 + x*C[l]>> 

6. DSolve[y’[x] == k y[x]/x, y[x], x] 

{{y[x] -> x k C[l]>> 

7. DSolvef{y*[x] == 1 + y[x] A 2, y[Pi/4] == Sqrt[3]}, 
y[x], x] 

{{y[x] -> Tan[Pi/12 + x]>> 

8. DSolve[ {n’[t] + lambda* n == 0, n[0] == n0>, n[t],t] 
{{n[t] -> n0 E A (- lambda*t>} 

9. DSolve[{R Q’[t] + Q[t]/C == 0, Q[0] == Qo},Q[t],t] 
{{Q[t] -> Qo/E-(t/(C*R))>> 

16.07 a) y = x (C + ln x) y = x (2 + ln x) 

b) y 2 + 2 x y - x 2 = K y 2 +2xy-x 2 =7 


c) K • x = (5 y 2 + x 2 ) 3 


21 x = (5 y 2 + x 2 ) 3 
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16. Differentialgleichungen erster Ordnung 


16.08* a) DSOLVEl_GEN(-y/x - 1, 1, x, y, c) 
y = x (LN (x) + c) 

DS0LVE1(-y/x - 1, 1, x, y, 1, 2) 

-LN(x) + y/x -2=0 
y = X LN (x) + 2 x 

b) DSOLVEl_GEN(y - x, x + y, x, y, c) 
y A 2 + 2 x y - x A 2 = 2 c 
DSOLVE1(y - x, x + y, x, y, 1, 2) 

-x A 2/2 + x y + (y A 2 -7 )/2 = 0 
y = V (2 x A 2 + 7) - x 
y = - V(2 x A 2 + 7) - X 

C) DS0LVE1_GEN(x A 2 - y A 2, 6 x y, x, y, c) 

(x A 2 + 5 y A 2) A 3 = k x 

DSOLVE1(x A 2 - y A 2, 6 X y, x, y, 1, 2) 

3 (x A 2 - 21 x A (1/3) + 5 y A 2) / (5 x A (l/3)) = 0 
y = V(5) V (x A (5/3) - 21) (—x) A (1/6) / 5 
y = - V(5) i/(x A (5/3) - 21) (-x) A (l/6) / 5 

16.09* a) DSolve[y'[x] == y[x]/x + 1, y[x], x] 

{{y[x] -> x C[l] + x Log[x]>} 

DSolve[{y'[x] == y(x]/x + 1, y[l] == 2}, y[x], x] 
{{y[x] -> 2 x + x Log[x]}} 

y[x_] = y[x]/.%[[!]]; Plot[y[x], {x, o, 1}]; 



DERIVE: siehe 16.08a 





16. Differentialgleichungen erster Ordnung 


b) gll = DSolve[(x + y[x]) y'[x] == x - y[x], y[x], x] 

-2 x - Sqrt[4 x 2 - 4 (-x 2 - 2 C[l])] 

{{y[x] ->-}, 

2 

-2 x + Sqrt[4 X 2 - 4 (-x 2 - 2 C[1])} 

<y[x] ->->} 

2 

DSolve[{(x + y[x]) y’[x] == x - y[x], y[l] == 2), 
y[x], x] 

-2 x + Sqrt[4 X 2 - 4 (-7 - X 2 )] 

{{y[x] ->->} 

2 

y[x_] = y[x]/.%[[l]]; Plot[y[x], {x, -2, 6>]; 


\ 

. 

\ 


mgm 


-2 -1 



\ / 



-2 

^ 9 ^ 4 6 



c) DSolve[6 x y[x] y’[x] == y[x] A 2 - x A 2, y[x], x] 
{{y[X] -> -(- X A 2/5 + X A (1/3) C[ 1] ) A (1/2)}, 

{y[X] -> (- X-2/5 + X A (1/3) C[l]) A (l/2)>> 

DSolve[{6 x y[x] y'[x] == y[x] A 2 - x A 2, y[l] == 2>, 
y[x], x] 

21 x 1 / 3 x 2 

{{y[x] -> Sqrt(-]}> 

5 5 

y[x_] - y[x]/.%[[!]]; Plot[y[x], {x,0,10}]j 
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16. Differentialgleichungen erster Ordnung 


16. 10 

a ) 

x 2 - y 2 = K 



b) 

x 3 - 3 x y + y 3 = 

c 

16. 11 

a) 

3 x 2 H 

H 4 x y 

- 4 

y 2 = C 

b) 

2 x 4 y 2 - 3 x 2 + 2 

y 3 = C 

16.12 

a) 

x 3 e y 

- 3 y = 

C 


b) 

x cos y = K 


16.13 

a) 

y = e' 

’ x3 [c + 

1 " 

{3 (x - 2) 

dx 

1 



b) 


(x 3 + 2 

x 2 

+ K x) 




16.14 

a) 

y = 4 

e _3x + 

x - 

2 

b) 

y - 5 - 2 

1 3 cos x 3 

cos 2 X 

16.15 

a) 


e 2x _ e 

X 

1 

b) 

y = K e" x + cos2x 

+ 2 sin2x 

y - * 



2 


5 

16.16 

a) 

y = c 

e sin x 

- 8 

(1 + sin x) 




b) 

y = k 

e~ X + 2 

sin x - cos 

X 



16.17 

a) 

y = K 




b) 

y = 1,12 x + x sin x 


sin x 


16.18 a) y = - cos x 


0,261080 





17. Differentialgleichungen zweiter Ordnung, 

die sich auf Differentialgleichungen erster Ordnung 
zurückführen lassen 

x 4 x 3 x 2 

17.01 a) y ■ } - § +| + C 1 x + C 2 

X = 0 0 = c 2 

X = 1 1 = 1/4 - 1/3 + 1/2 + C x =» C 1 = 7/12 



MATHEMATICA: 

DSolve[y"[x] == 3 x A 2 - 2 x + 1, y[x], x] 

{{y[x] -> X A 2/2 - X A 3/3 + X A 4/4 + C[l] + X C[2]>> 
DSolve[{y” [x] == 3 x A 2 - 2 x + 1, y[0] == 0, y[l] == 1}, 

y[x], *] 

{{y[X] -> (7 X)/12 + X A 2/2 - X A 3/3 + X A 4/4>} 

Plot[7 x/12 + x A 2/2 - x A 3/3 - x A 4/4, {x, -5, 5}, 
PlotRange->{-6, 1}]; 



x = 0 0 = c 2 

X = 1 1 = 2 sin 1 - COS 1 + C 2 =» C ± = - 0,14264 

y = 2 sin x - x cos x - 0,14264 x 

MATHEMATICA: 

DSolve[y"[x] = x Cos[x], y[x], x] 

{{y[x] -> C[l] + x C[2] - x Cos[x] + 2 Sin[x]>} 
DSolve[{y’*[x] == x Cos[x], y[0] == 0, y[l] == 1}, 
y[x], x] 

{{y[x] -> -0.14264 x - 1. x Cos[x] + 2. Sin[x]>> 
Plot[-0.14264 x - x Cos[x] + 2 Sin[x], {x, -5, 5}]; 
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17.03 a) vT y 
2 di 


S-p 

'!• 


Probe: y ' 


dp _ dx 
P 2 


dx + O 2 

, 1/X 


p"-iI 

ln p = - j + ln 0 
singuläre Lösung: y < 


DERIVE und MATHEMATICA liefern keine Lösung, weil sie 


J‘ 


dx nicht direkt lösen können. 


b) y = Cj^ln x + C 2 singuläre Lösung: y = k 

MATHEMATICA: 

DSolve[x y”[x] + y* [x] == o, y[x], x] 

{{y[x] -> C[i] + C[2] Log[x]>> 

DSolve[{x y"[x] + y’[x] == 0, 

y[1] == 2, y* [l] == l}, y[x], X] 

{{y[x] -> 2 + Log[x]}} 

Plot[2 + Log[x],{x,-2,4}]; 
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b) y' = p(y) 


y # ' = Pdp/dy = p.p' 
singuläre Lösung: y = 


! ~2 ’ 


MATHEMATICA: 

DSolve[y[x] y'*[x] == y'[x] A 2, y[x], x] 

DSolve::dnim: 

Built-in procedures cannot solve this differential 
equation. 

DSolve[y[x] y"[x] == y'[x] A 2, y[x], x] 

"Ersetzung y' = p, y " = p-dp/dy = p.p' 

(y ist unabh. Var.!)"; 

DSolve[y p[y] p'[y] == p[y] A 2, p[y], y] 

{(P[y] -> E a (C[ 1] + Log[y])}) 

"p = y' wieder einsetzen"; 

DSolve[y»[x] == Exp[Log[c y[x]]], y[x], x] 

{{y[x] -> E A (c x) C[l]}} 

Plot[2/E A 2 Exp[2 x], {x, -3, 1}, PlotRange->{0, 1> J; 



-3 


-2 


-1 


0 


1 



18. Differentialgleichungen zweiter Ordnung 
mit konstanten Koeffizienten 

18.01 a) y = C ± e~ x + C 2 e _2x b) y = 7 e 3x - 5 e 4x 

18.02 a) y = e X ^ x + C 2 j b) y = e _3x x + 2 J 

18.03 a) y = e x ^ cos 2x + C 2 sin 2x| 
b) y = e -2x |cos x + 2sin xj 


18.04 a) y = e 

b) y = 8,83340 10 _ ‘*e J ’ x sin 2 x 


cos 1,936 x + C 2 sin 1,936 xj 
4 3-x , 


18.05 (18.01 a) mit MATHEMATICA) 

DSolve[y*’[x] + 3 y'[x] + 2 y[x] == 0, y[x], x] 


{{y[x] -> 


,2 x 


->> 


18.06 (18.02 a) mit MATHEMATICA) 

DSolve[y’’[x] - 2 y'[x] + y[x] == 0, y[x], x] 

{{y[x] -> E X C[l] + E X x C[2]}> 

18.07* (18.03 a) mit MATHEMATICA) 

DSolve[y*'[x] - 2 y’[x] + 5 y[x] == 0, y[x], x] 
{{y[x] -> E x C[2] Cos[2 x] - E x C[l] Sin[2 x]}> 

18.08* (18.04 a) mit MATHEMATICA) 

DSolve[y'*[x] + 5 y»[x] + 10 y[x] == 0, y[x], x] 
{{y(x] -> E^- 5 - 1 x >/ 2 C [l] 

+ e ((-5 + I Sqrt(15]) x)/2 c[2]}} 


18.09 (18.01 b) mit DERIVE) 

DSOLVE2_IV(-7,12,0, x, 0,2,1) 





18. Lineare Differentialgleichungen zweiter Ordnung 


18.10* (18.02 b) mit DERIVE) 

DSOLVE2_IV(6, 9, 0, X, 0, 2, 0) 
e -3 X (6 x + 2) 


18.11* (18.03 b) mit DERIVE) 

DS0LVE2_IV( 4 , 5, 0, X, 0, 1, 0) 
e _2x (COS(x) + 2 SIN(x)) 


18.12* (18.04 b) mit DERIVE) 

DSOLVE2_BV(-6, 13, 0, X, 0, 0, 3, 
8.8334010 -4 -e 3 ' X SIN(2 x) 



2x 


18.13 a) y h = cosv'j? x + C 2 sim/fs x 

y = ^ cosvl) x + C 0 sinv'j? x - 

b) y h = C x e 

y = c i e 

18.14 a) y = e 


C 0 e‘ 

C. e 


2x 3 2 1 , 


X - 2 X + 4 

1/42X A „ ö -8,42x y _ b . e : 

y = - 1,25 e x + C x e 1,42x + C 2 e ~ 8,42x 
b ) y h = e X/2 (c x cos ^ X + c 2 sin ^ x] y p = b e 3x 

y = e x/2 [c n cos 1,66 x + C, sin 1,66 x] + 0,222 e 3: 






Lineare Differentialgleichungen zweiter Ordnung 


18.15 a) y = r + C_ 


a x + bx + cx + dx + < 


y = O + c 2 e -------- 

| e^^ bleibt unbestimmt, + e^^ = C j 

b) y h = c ± + C 2 e" 5x y p = a x 3 + b x 2 + c x 

y = C + C. 

f d bleibt unbestimmt, + d = C 


-5x , x_ 13x 113x 

15 50 “ 125 


18.16 a) y h ■ 

y ; 

b) y h ■ 

y ; 

18.17 a) y ■ 


l C l + C 2 X J 

(cj + C 2 X + 2 X 2 ) 

( C 1 + C 2 *) 

K + C 2 X + 5 X 1 

(c. cos/? x + C, sin/? xl 


y = e _2x |c x cos/? x + C 2 sin/? xj + j x ■ 


121 

“49 


b) y h ■ 


+ G x cos x + H cos x 

y - «' 1 ( » «/H . , 

- 4 x sin x + 3 x cos x - 8 cos x + 5 sin x 
h = C e 5x + C e* 4x y = e x (a x + b) + c 


= C. e 5x + C. 


-4x xe e 3 

20 400 20 


b) y t 
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y h = e cos 3x + C 2 sin 3xJ 

Yp = a sin 2x + b cos 2x 

y = e -x (c, cos 3x + c o sin 3x] + sin 2x “ t? cos 2x 


13 


C, cos 4x + C_ sin 4x 


sin 4x - == cos 4x + öc cos 3x ~ Tk sin 3x 


!8.21 Y h = c ± e 
y = C, e 3: 


C 2 e 
C. e" 


y = b x e 


18.22 y h = e + C 2 e y p = b x e 

y = - 1 e 4x - M e" 2x -Ix e _2x 

y 9 e 9 e 6 Xe 

18.23 y h = C 1 cos 3x + C 2 sin 3x 

Yp = A x cos x + B x sin x + C cos x + D sin x 
y = 9 ^ sin 3x + I x cos x - II sin x 

18.24 y h = C 1 e 5x + C 2 e _x 

Man rechnet mit besserer Übersicht, wenn man 
Yp in Yp = Yp + Yp aufspaltet und jeden Teil 
für sich allein rechnet. 

Yp = Ax 2 +Bx + C 

y = c e 5x + c 


-x l 2 


13 


x cos x - 


169 


5 X + 

29 

338 


125 26 

cos x 


18.25 (18.19 mit MATHEMATICA) 

DSolve[y* *[x] + 2 y*[x] + 10 y[x] == 0,y[x],x] 
C[2] Cos[3 x] C[l] Sin[3 x] 

{{y[ x ] ->- z - n- }> 
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18. Lineare Differenentialgleichung zweiter Ordnung 


11 = DSolve[y»'[x] + 2 y* [x] + 10 y[x] == 3 Sin[2 x], y[x],x] 
C[2] Cos[ 3 X] 


(— + —) 

E” 52 52 

C[l] Sin[3 X] 

((-6 + 6 1) Cos [2 x] + (9 - 9 I) Sin [2 x])---}} 


Nun hilft uns der Befehl ExpandAll weiter. 

ExpandAl1 [ausdruck] multipliziert Produkte und Potenzen mit 
positiven ganzen Hochzahlen in jedem Teil von ausdruck aus. 

ExpandAll [ausdruck, Trig -> True] behandelt trigonometrische 
Funktionen wie rationale Funktionen mit Exponentialaus- 
drücken und entwickelt sie entsprechend. 


12 = ExpandAll[11[[1]]] 

-3 Cos[2 x] C[2] Cos[3 x] 9 Sin[2 x] C[l] Sin[3 x] 

{y[x] ->-+---+---> 

13 E x 26 E x 

Nach dieser Einleitung behandeln wir unser eigentliches Problem, 
das Lösen der Differentialgleichungen mit Randbedingungen. 

13 = DSolve[{y*'[x] + 2 y*[x] + 10 y[x] == 3 Sin[2 x], 
y[0] == 0, y[Pi/2] == 3/13 >, y[x], x] 

3 Cos[3 x] 


Hy[x] -> 


13 E 


+ — + — (-6 + 6 I) Cos[2 x] + (9 - 9 I) Sin[2 x] }} 

'■52 52' 1 > 

14 = ExpandAl1[13[[1]]] 

—3 Cos[2 x] 3 Cos[3 X] 9 Sin[2 X] 

-+-j— +- 

13 E 26 


{y[x] -> • 
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18.26* (18.21 mit MATHEMATICA) 

11 = DSolve[{y’’ [x] + 4 y’[x] - 21 y[x] == 2 Exp[3 x], 

y[1.5] == l, y*[1.5] == -2>, y(x], x] 

83538.1 3 x E 3 X X 

{ (Y[X] -> - 0.314446 E +->} 

E' X 5 

y[x_] = y[x]/.ll[[l]]; 

Flot[y[x], {x, 1, 2}, PlotRange -> {0, 5> 



1 1.2 1.4 1.6 1.8 2 


18.27* (18.23 mit MATHEMATICA) 

11 = DSolve[(y" [x] + 9 y[x] == x Cos[x] - 3 Sin[x], 
y[0 ] == 0 , y»[0] == 0}, y[x], x] 

{{y[x] -> (Cos[2 x] + I Sin[2 x]) 


-51 x (Cos[3 x] - I Sin[3 x]) 

•(-(Cos [3 x] - I Sin[3 x)) +-) 


48 24 

+ _ 

Es folgen noch Ausdrücke, die mehr als eine halbe A4-Seite 
Platz brauchen. Daher: 12 = ExpandAll[11[[1]]] 

{y[x] -> (x Cos[2 x] Cos[3 x])/12 + (x Cos[3 x] Cos[4 x])/24 
+ (5 Cos[3 x] Sin[2 x])/24 + (7 Sin[3 x])/96 - 

- (5 Cos[2 X] Sin[3 x])/24 + (13 Cos[4 X] Sin[3 X] )/96 

+ (x Sin[2 x] Sin[3 x])/12 - (13 Cos[3 x] Sin[4 x])/96 

+ (x sin[3 x] sin[4 x])/24} 

y[_] = y[x]/.13[[l]]; Plot[y[x], {x, -2 Pi, 2Pi}]; 
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18.30* (18.24 mit DERIVE) 

DSOLVE2(“4, -5, X A 2 - 2 + X SIN(x), X, Cl, c2) 

Cl #e 5x + c2 #e” X + 1/338 (26 x - 29) COS(x) 

-1/338 (39 x + 2) SIN (x) - 0.008 (25 x - 40 x - 8) 2 
DSOLVE2_IV(-4, -5, X A 2 - 2 + X SIN(x), x, 0, 0.75, 1.2) 
0.263465 #e 5 X + 0.508333 #e" X 
+ 0.00295857 (26 X - 29) COS(x) 

- 0.00295857 (39 X + 2) SIN(x) 

- 0.008 (25 X 2 - 40 X - 8 ) 



18.31 (1. Fall (Seite 220) mit MATHEMATICA) 

Remove["Global**"] 

dg = DSolve[{9 y* *[t] + 6 y*[t] ♦ 82 y[t] ” 0» 
y[0] == 0, y'[0] == 3>, y[t], t] 
{{y[t] -> Sin[3 t]/E-(t/3)}> 
y = y[t]/.dg[[i]]; 

Plot[{y, Exp[-t/3], -Exp[-t/3]}, (t, -1, 2Pi}]; 







Lineare Differentialgleichungen zweiter Ordnung 



FindMinifflum[-y, {t, 0}] {-0.845039, {t -> 0.486713}} 

FindNinimum[y, {t, 1.5}] {-0.596045, {t -> 1.53391}} 

k = D[y, t] (3*Cos[3*t])/E A (t/3) - Sin[3*t]/(3*E A (t/3)) 

k /.t-> O 3 

t := 0.8; 

Solve[z - y == k (x - t), z] {{z -> 2.01082 - 1.86682 x}} 


18.32 (2. Fall (Seite 220) mit MATHEMATICA) 

Remove["Global**"] 

dg = DSolve[ {y* *[ t ] + 3 y’[t] + 2 y[t] == 0, 
y[°] == i, y’[0] == -1), y[t], t] 
{{y[t] -> e a (- t)}} 
y = y[t] /. dg[[ 1 ]]; 

Plot[y, {t, 0, 2Pi}]; 
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18.33 (3. Fall (Seite 221) mit MATHEMATICA) 

Remove["Global'*"] 

dg = DSolve[ {y" [t] + 4 y' [t] + 4 y[t] == 0, 
y[°] == i# y’[0] == -4}, y[t], t] 
{{y[t] -> E-(-2*t) - (2*t)/E-(2*t)}} 
y = y[t]/.dg([l]]; 

Plot[y,{t, 0, 2Pi)]; 


0.2 

0.15 

0.1 

0.05 

- 0.05 

- 0.1 



4 5 6 


FindMinimum[y, {x, 0.8}] 

{-0.135335, {X -> 1}} 

k = D[y,t] 

-4/E* (2 t) + (4 t) /E A (2 t) 
t := 0; 

Solve[z - y == k (x - t), z] 

{{Z -> 1 - 4 x}} 
t := 0.8; 

Solve[z - y == k (x - t), z] 

{{Z -> 0.00807586 - 0.161517 X}} 

Minimum ( 1 | -0,1353 ) 

Gleichung der Tangente in ( 0| 1) : _ y = 1 - 4 x 


Gleichung der Tangente in( 0,8 | y ) 


y = 0,0081 - 0,1615 x 




19. Anwendung der Vektorrechnung 

19.01 

1 : x := 1 + lambda 
2 : y := 2 lambda 

3: z := 2 + 3 lambda 

4: (x + 2) A 2 + (y - 3) A 2 + (z - 5) A 2 = 36 
5: 14 lambda A 2 - 24 lambda + 27 = 36 
6 : lambda = -0.316548 

8 : lambda = 2.03083 

9: "Schnittpunkt Sl:" 

10: Xl := 0.683452 

11: yl := -0.633096 

12: zl := 1.05035 

13: "Schnittpunkt S2:" 

14: x2 := 3.03083 

15: y2 := 4.06166 

16: z2 := 8.09249 

17: "Länge der Sehne S1S2:" 

18: d := ✓((xl " X2) A 2 + (yl - y2) A 2 + (zl - z2) A 2) 8.78308 

20: "Tangentialebenen:" 

21: m := [-2, 3, 5] 

22: sl := [xl, yl, Zl] [0.683452, -0.633096, 1.05035] 

24: s2 := [x2, y2, z2] [3.03083, 4.06166, 8.09249] 

26: xv := [x, y, z] 

27: "Ebene durch Sl:" 

28: (xv - sl)-(m - sl) = 0 

29: -2.68345 X + 3.63309 y + 3.94965 Z - 0.0144056 = 0 

30: "Ebene durch S2:" 

31: (xv - s2 ) •(m - s2) = 0 

32: -5.03083 x - 1.06166 y - 3.09249 z + 44.5856 = 0 

19.02 

1: "Ebene E2:" 

2: 10 + 2 5 + 1 = c c = 21 

4: "E2: x + 2 y + z = 21" 

5: "Abstand der beiden Ebenen = Abstand PE1 = 2r:" 

6 : 2 r = (10 + 2 5 + 1 - 8 ) / ^(1-2 + 2 A 2 + 1 A 1) 

7: r = 13 ✓(6)/12 r = 


2.65361 
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9: "Ebene E3 (parallel E1,E2) durch den Kugelmittelpunkt:" 

10 : "x + 2 y + z = e" 

11: "P hat von E3 den Abstand r:" 

12: (10 +2 5 + 1 -e) / /(1*2 + 2*2 + 1*2) = 2.65361 

13: e = 14.5 

14: "E3: x + 2 y + Z = 14.5" 

15: "Mittelpunkt M(ml,m2,m3) ist Schnitt von E3 mit g" 

16: "g: xv = (5, 3, 2) + lambda*(-l, -1, -7)" 

17: (5 - lambda) + 2 (3 - lambda) + (2 - 7 lambda) = 14.5 

18: lambda := -0.15 

20: m := [5, 3, 2] + lambda [-1, -1, -7] [5.15, 3.15, 3.05] 

22: "K(5.15, 3.15, 3.05; 2.65361)" 

23: (x - 5.15)*2 + (y - 3.15)*2 + (z - 3.05)*2 = (2.65361)*2 

24: "Kontrolle: Abstand ME1:" 

25: (5.15 + 2 3.15 + 3.05 - 8) / /(1*2 + 2*2 + 1*2) 

26: 13 /(6)/12 2.65361 

19.03 

1: "Koordinaten von TI:" 

2: Solve(2 5 + 3 y + 6 2 = 24, y) y= 2/3 

4: tl := [5, 2/3, 2] 

5: "Normalvektor der Ebene:" 

6 : n := [2, 3, 6] 

7: "Gerade durch Tl und M:" 

8 : tl + lambda n 

9: [2 lambda + 5, 3 lambda + 2/3, 6 lambda + 2] 

10: m := [2 lambda + 5, 3 lambda + 2/3, 6 lambda + 2] 

11: pl := [8, -5, -4] 

12: "Tl auf K:" 

13: (tl - m)*2 = r*2 

14: 49 lambda*2 = r*2 

15: "Pl auf K:" 

16: (pl - m)*2 = r*2 

17: (pl - m)*2 = 49 lambda*2 

18: (441 lambda*2 + 846 lambda + 694)/9 = 49 lambda*2 

19: Solve(#18, lambda) lambda = -347/423 

21: Solve(49 lambda*2 = r*2, r); 

22: r = 7 lambda r = 5.74231 
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23: m 

24: [2 lambda + 5, 3 lambda + 2/3, 6 lambda + 2] 

25: [3.35933, -1.79432, -2.92198] 

26: "Kugel K(3.35933, -1.79432, -2.92198; 5.74231)" 

27: "Tangentialebene im Punkt PI:" 

28: xv := [x, y, z] 

28: (xv - pl)(m - pl) = 0 

29: (3 x (2 lambda - 3) + y (9 lambda + 17) + 18 z (lambda + 1) 

+ 69 lambda + 229)/3 = 0 

30: -4.64066 x + 3.20567 y + 1.07801 Z + 57.4657=0 

31: "Hinkel zwischen den beiden Ebenen:" 

32: nl := [2,3,6] 

33: n2 := [-4.64066, 3.20567, 1.07801] 

34: AC0S(nl n2 / /(nl A 2 n2 A 2))/deg 80.2549 

19.04 

1: m := [3, -7, -5] + lambda [1,2,3] 

2: pl := [10, 3, 1] 

3: p2 := [4, -5, 8.4833] 

4: "Pl auf K:" 

5: (pl -m) A 2 = r A 2 

6 : "P2 auf K:" 

7: (p2 - m) A 2 = r A 2 

8 : (pl - m ) A 2 = <p 2 - m ) A 2 

9: 14 lambda A 2 - 90 lambda + 185 = 10 A (-8) (1.4 10 A 9 lambda A 2 

-9.08998 10 A 9 lambda + 1.86799*10 A 10) 

10: Solve(#9, lambda) lambda = 1.99975 

12 : lambda := 2 

13: m; [5, -3, 1] 

15: (pl - m) A 2 = r A 2 

16: Solve(#15, r) r = 7.81024 

18: "Kugel K(5, -3, 1; 7.81024)" 

19: xv := [x, y, z] 

20: "Tangentialebene in Pl:" 

21: (xv - pl)•(m- pl) =0 -5 x - 6 y + 68 = 0 

23: "Tangentialebene in P2:" 

24: (xv - p2)-(m - p2) =0 x + 2 y - 7.4833 Z + 69.4830 = 0 

26: "Hinkel zwischen den beiden Ebenen:" 

27: nl := [-5, - 6 , 0] 

28: n2 := [1, 2, -7.4833] 

29: ACOS(nl n2 / /(nl A 2 n2 A 2)) / deg; 


106.181 




21. Simplexmethode 


Das "händische" Lösen von Aufgaben der linearen Optimierung ist 
nur für ganz kleine Aufgaben sinnvoll. 

Wir. bringen daher ein einfaches QBASIC-Programm für die Berechnung 
von nichtentarteten Lösungen, bei dem die einzelnen Schritte 
studiert werden können. 

Für Aufgaben der Praxis können Sie die Ausgabe der Zwischener¬ 
gebnisse abstellen.. 


DECLARE SÜB druck (a(), n!, m!) 

' Simplex-Methode 

CLS 

OPEN "simp.dat" FOR OUTPUT AS #1 
DIM a(10, 20), b(10, 20), x(10) 

READ n, m 

FOR i = 1 TO n + 1 

FOR j = 1 TO n + m + 1 : READ a(i, j) : NEXT j 
NEXT i 

CALL druck(a(), n + 1, n + m + 1) 

FOR k = 1 TO n 

s = 1: gr = a(n + 1, 1) 

FOR j = 2 TO n + m 

IF a(n + 1, j) > gr THEN 

gr = a(n +1, j): s = j 'Pivotspalte 

END IF 
NEXT j 

IF gr <= 0 THEN EXIT FOR 

FOR i = 1 TO n: a(i, n + m + 2) = 1E30: NEXT i 
FOR i = 1 TO n + 1 

IF a(i, s) > 0 THEN a(i, n + m + 2) = a(i, n + m + l)/a(i, s) 
NEXT i 

z = 1: kl = a(l, n + m + 2) 

FOR i = 2 TO n 

IF a(i, n + m + 2) < kl THEN 
kl = a(i, n + m + 2): z = i 
END IF 

IF kl = 1E30 THEN PRINT "keine Lösung” : END 
NEXT i: PRINT #1, z, s 'Pivotzeile, Pivotspalte 

' Pivotzeile 

FOR j = 1 TO n + m + 1: b(z, j) = a(z, j) / a(z, s): NEXT j 
' Elimination 

FOR i = l TO n + 1: f= a(i, s) 

IF i <> z THEN 

FOR j=lTOn+m+l 

b(i, j) = a(i, j) - f * b(z, j) 

NEXT j 
END IF 
NEXT i 



94 


21. Simplexmethode 


FOR i = 1 TO n + 1 

FOR j =lTOn+m+l 
a(i, j) = b(i, j) 

NEXT j 
NEXT i 

CALL druck(a () , n+1, n+m+2) 

NEXT k 

' Lösungen 

FOR j = 1 TO n + m 

IF a(n + 1, j) <0 THEN 
x(j) = 0 : p = p +1 
ELSE 

FOR i = 1 TO n: IF a(i, j) = 1 THEN x(j) = a(i, n + m + 1) 

IF x(j) < 0 THEN PRINT "keine Lösung" : END 
NEXT i 
END IF 
NEXT j 

If p < m THEN PRINT "entartete Lösung! " : END 
FOR j = 1 TO m 

PRINT #1, "x("; j; ") = "; USING "####.###"; x(j) 

NEXT j 

PRINT #1, : PRINT #1, "Z roax = "; -b(n +1, n + m + 1): PRINT #1, 

SUB druck (a(), n, m) 

FOR i = 1 TO n 
FOR j = 1 TO m 

PRINT #1, USING "####.##"; a(i, j); 

NEXT j: PRINT #1, 

NEXT i 
PRINT #1, 

END SUB 

DATA 3,4, 2,1,0,1,1,0,0,10 

DATA 1,-1,7,-1,0,1,0,30 
DATA 1,1,3,0,0,0,1,25 
DATA 5,6,1,2,0,0,1,0 

2.00 1.00 0.00 1.00 1.00 0.00 0.00 10.00 

1.00 -1.00 7.00 -1.00 0.00 1.00 0.00 30.00 

1.00 1.00 3.00 0.00 0.00 0.00 1.00 25.00 

5.00 6.00 1.00 2.00 0.00 0.00 1.00 0.00 

1 2 

2.00 1.00 0.00 1.00 1.00 0.00 0.00 10.00 10.00 

3.00 0.00 7.00 0.00 1.00 1.00 0.00 40.00 100.00 

-1.00 0.00 3.00 -1.00 -1.00 0.00 1.00 15.00 25.00 

-7.00 0.00 1.00 -4.00 -6.00 0.00 1.00 -60.00 0.00 

3 3 

2.00 1.00 0.00 1.00 1.00 0.00 0.00 10.00 100.00 

5.33 0.00 0.00 2.33 3.33 1.00 -2.33 5.00 5.71 

-0.33 0.00 1.00 -0.33 -0.33 0.00 0.33 5.00 5.00 

-6.67 0.00 0.00 -3.67 -5.67 0.00 0.67 -65.00 -60.00 

X( 1 ) = 0.000 

X( 2 ) = 10.000 

X( 3 ) = 5.000 

X( 4 ) = 0.000 


Z max 


65 
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21.01 







2.00 

3.00 10.00 

1.00 

0.00 

0.00 

300.00 


1.00 

- 2.00 1.00 

0.00 

1.00 

0.00 

30.00 


- 5.00 

8.00 3.00 

0.00 

0.00 

1.00 

100.00 


3.00 

4.00 10.00 

0.00 

0.00 

0.00 

0.00 


1 

3 






0.20 

0.30 1.00 

0.10 

0.00 

0.00 

30.00 

30.00 

0.80 

- 2.30 0.00 

- 0.10 

1.00 

0.00 

0.00 

30.00 

- 5.60 

7.10 0.00 

- 0.30 

0.00 

1.00 

10.00 

33.33 

1.00 

1.00 0.00 

- 1.00 

0.00 

0 . 00 - 

- 300.00 

0.00 

2 

1 






0.00 

0.88 1.00 

0.13 

- 0.25 

0.00 

30.00 

150.00 

1.00 

- 2.88 0.00 

- 0.13 

1.25 

0.00 

0.00 

0.00 

0.00 

- 9.00 0.00 

- 1.00 

7.00 

1.00 

10.00 

100.00 

0.00 

3.88 0.00 

- 0.88 

- 1.25 

0 . 00 - 

- 300 . 00 - 

• 300.00 

1 

2 






0.00 

1.00 1.14 

0.14 

- 0.29 

0.00 

34.29 

34.29 

1.00 

0.00 3.29 

0.29 

0.43 

0.00 

98.57 

100.00 

0.00 

0.00 10.29 

0.29 

4.43 

1.00 

318.57 

100.00 

0.00 

0.00 - 4.43 

- 1.43 

- 0.14 

0 . 00 - 

- 432.86 

- 77.42 

x( 1 ) = 

98.571 






X( 2 ) = 

34.286 






X( 3 ) = 

0.000 






Z max = 

432.8571 






21.02 







2.00 

3.00 10.00 

1.00 

0.00 

0.00 

300.00 


1.00 

- 2.00 1.00 

0.00 

1.00 

0.00 

30.00 


- 5.00 

8.00 3.00 

0.00 

0.00 

1.00 

100.00 


1.00 

- 4.00 15.00 

0.00 

0.00 

0.00 

0.00 


1 

3 






0.20 

0.30 1.00 

0.10 

0.00 

0.00 

30.00 

30.00 

0.80 

- 2.30 0.00 

- 0.10 

1.00 

0.00 

0.00 

30.00 

- 5.60 

7.10 0.00 

- 0.30 

0.00 

1.00 

10.00 

33.33 

- 2.00 

- 8.50 0.00 

- 1.50 

0.00 

0 . 00 - 

- 450.00 

0.00 

X( 1 ) = 

0.000 






X( 2 ) = 

0.000 






X( 3 ) = 

30.000 






Z max = 

450 






21.03 







ConstrainedMax[ 30 xl + 

24 X 2 - 

15 X 3 , 





{ 2 x 1 - 

3 X 2 - 

10 x 3 <= 

30 , 




Xl + 

2 x 2 + 

8 x 3 <= 

30 , 




5 x 1 + 

18 x 2 - 

5 x 3 <= 

100 }, 




{Xl,x 2 

!,X 3 > ]//N 




{ 608 . 891 , 

{Xl -> 20 . 7839 , X 2 

-> 0 . 0956023 , 

X 3 -> 1 

. 12811 }} 
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21.04 

ConstrainedMax[ 8x1 - x2 + 30x3, 

{ 2 X 1 - 1 X 2 + 2 X 3 <= 80, 

3x1 + 1 x2 + 2x3 <=120 >, 
{xl,x2,x3> ]//N 

{1480., {Xl -> 0, x2 -> 20., X 3 -> 50.}} 


21.05 

ConstralnedMax[50xl + 70x2 + 20x3, 

{ 2x1 - x2 + 12x3 < 

3x1 - 4 x2 + 20x3 < 
xl x3 < 

xl + 2 x2 < 

5x1 - 6X3 < 

2x1 + 3x3 < 

{xl,x2,x3> ]//N 

{423.2, {Xl -> 3.68, x2 -> 3.16, x3 -> 0.9}} 


21.06 

ConstrainedMax[150x1 + 70x2 + 40x3, 

{ 2x1 - x2 + 12X3 <= 150, 

3x1 - 4 x2 + 20x3 <= 250, 

Xl - x3 <= 30, 

Xl + 2 X2 <= 100, 

5x1 - 6x3 <= 130, 

2x1 + 3x3 <= 120 }, 

{xl,x2,x3} ]//N 

{8092., {Xl -> 36.8, x2 -> 31.6, x3 -> 9.}} 

21.07 

ConstrainedMin[-2xl + x2 - 10x3, 

{ 4x1 - 2X2 + 2x3 <= 3, 

Xl + 5x2 + 15x3 <= 65 }, {Xl,x2,x3} ]//N 

{-34.125, {Xl -> 0, x2 -> 2.125, x3 -> 3.625}} 

21.08 

ConstrainedMax[-2xl + x2 + 10x3, { 4x1 - 2x2 + 2x3 <= 3, 

xl + 5x2 + 5x3 <=65 }, {xl,x2,x3} ]//N 

{78.25, {Xl -> 0, X2 -> 5.75, x3 -> 7.25}} 

21.09 

ConstrainedMin[ 2x1 + x2 - 10x3, 

{ 4X1 - 2x2 + 2x3 <= 3, 

3X1 + 5X2 + 15X3 <= 65 }, {Xl,x2,x3} ]//N 

{-34.125, {Xl -> 0, x2 -> 2.125, X3 -> 3.625}} 

21.10 

ConstrainedMax[ 5x1 + 6x2 - 10x3, 

{ 4xl - 2x2 + 2x3 <= 3, 

Xl + 5x2 + 15x3 <= 65 }, 

{xl,x2,x3} ]//N 

{103.045, {Xl -> 6.59091, x2 -> 11.6818, x3 -> 0}} 


= 15, 

= 25 , 

= 3, 

= 10 , 

= 13, 

= 12 } 
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22. Netzplantechnik 


Die beiden Programme sind den Im Verlag Oldenbourg erschienenen 
Büchern 

Blaha / Schärf; "Programmiertraining BASIC", Band 2, 
bzw. 

Blaha / Schärf / Haselberger; "Programmiertraining PASCAL", 
Band 2 

entnommen. 


DEFINT A, E, I-J, V ' Programm CPM 

' Daten: " ", Vorgangsanfang,-ende,-dauer. Abschließen mit * 

' Reihung: vorhergehende Vorgänge müssen beendet sein 

CLS 

DIM Vor$(50), i(50), j(50), d(50), FAZ(50), FEZ(50) 

DIH SAZ(50), SEZ(50), GP(50), FZ(30), SZ(30) 
v$ = "### 

a$ = "Anf- End- Vorgs- früh. spät. früh. spät. Gesamt-" 
b$ = " Ereignis dauer Anfangszeit Endzeit Pufferzeit" 
e$ = " ## ### ###" 

c$ = "Ereignis FZ SZ" 

v = 0: AnzEr = 0 'Einlesen der Vorgänge 

DO 

READ Vor$: IF Vor$ = "*" THEN AnzVg = v: EXIT DO 
v = v + 1: Vor$(v) = Vor$: READ i(v), j(v), d(v) 

IF j(v) > AnzEr THEN AnzEr = j(v) 

LOOP 


'Reihung der Vorgänge 

FOR v = 1 TO AnzVg - 1: FOR w = v + 1 TO AnzVg 

IF j(w) = i(v) THEN PRINT "Falsche Reihung der Vorgänge !": END 
NEXT w, v 


'Vorwärtsrechnung 

FOR e = 1 TO AnzEr: FZ(e) = 0: NEXT e: F = 0 
FOR v = 1 TO AnzVg 

FAZ(v) = FZ(i(v)): FEZ(v) = FAZ(v) + d(v) 

IF FEZ(v) > FZ(j(v)) THEN FZ(j(v)) = FEZ(v) 

IF FEZ(v) > F THEN F = FEZ(v) 

NEXT v 


'Rückwärtsrechnung 

FOR e = 1 TO AnzEr: SZ(e) = F: NEXT e 
FOR v = AnzVg TO 1 STEP -1 

SEZ(v) = SZ(j(v)): SAZ(v) = SEZ(v) - d(v) 

IF SAZ(v) < SZ(i(v)) THEN SZ(i(v)) = SAZ(v) 

NEXT V 


'Tabelle der Vorgänge 

CLS : PRINT a$: PRINT b$ 

FOR v = 1 TO AnzVg 

GP(v) = SAZ(v) - FAZ(v) 

PRINT USING v$; i(v); j(v); d(v); FAZ(v); FEZ(v); 

PRINT USING "### "; SAZ(v); SEZ(v); GP(v); 

IF GP(v) THEN PRINT ELSE PRINT " K" 

NEXT V 
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' Tabelle der Ereignisse 

PRINT : PRINT TAB(10); c$; 

FOR e = 1 TO AnzEr 

PRINT TAB(10); USING e$; e; FZ(e); SZ(e); 

NEXT e 


DATA " ",1,2,4, 
DATA " ",2,7,5, 
DATA " ",5,6,7, 
DATA " ",7,8,6, 


",2,3,1, 

”,3,4,3, 

",5,7,6, 

",8,9,8, 


" ",2,5,2 
" ",4,8,2 
" ",6,7,1 
«*« 


Wir erhalten die Ausgabe: 


Anf- End- 
Ereignis 
1 2 

2 3 

2 5 

2 7 

3 4 

4 8 

5 6 

5 7 

6 7 

7 8 

8 9 


Vorgs- früh, spät, 
dauer Anfangszeit 

4 0 4 

14 5 

2 4 6 

5 4 9 

3 5 8 

2 8 10 

7 6 13 

6 6 12 

1 13 14 

6 14 20 

8 20 28 


Ereignis 

1 

2 

3 

4 

5 

6 

7 

8 
9 


FZ SZ 

0 0 

4 4 

5 15 

8 18 

6 6 

13 13 

14 14 

20 20 

28 28 


früh. spät. 
Endzeit 
0 4 

14 15 

4 6 

9 14 

15 18 

18 20 

6 13 

8 14 

13 14 

14 20 

20 28 


Gesamt- 
Pufferzeit 
0 K 

10 

0 K 

5 

10 

10 

0 K 

2 

0 K 

0 K 

0 K 


Beachten Sie: 

Für das Lösungsheft wird die Computer-Ausgabe 1 : V2 verkleinert. 
Um die Seitenbreite von 66 Buchstaben zu erhalten, haben wir in 
der Tabelle die Namen der Vorgänge weggelassen. 
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Program CPM; {PASCAL} 

uses crt; 

var i,j,D:array[l..30] of integer; 

{Anfangs-,Endereignis,Dauer d.Vorgangs} 
AnzVg,AnzEr,FE,e,v,w:integer; 

FAZ,FEZ,SAZ,SEZ,FZ,SZ: array[1..30] of integer; 
f: text; FN: string[20]; 

begin 

clrscr; 

write('Eingabedatei = ');readln (FN); writeln; 
assign (f,FN);reset(f); 

{Einlesen der Vorgänge} 

AnzEr:=0; v:=l; 
repeat 

readln(f,N[v],i[v],j[v],D[v]); 
if j[v]>AnzEr then AnzEr:=j[v]; 
v:=v+l; 

if eoln(f) then readln(f); 
until eoln(f); 

AnzVg:=v-l; 


for v:=l to AnzVg-1 do 
for w:=v+l to AnzVg do 
if j[w]=i[v] then begin 

writeln(' Falsche Reihung !'); 


{Reihung} 


readln; exit; end; 


{Vorwärtsrechnung} 

FE:=0; 

for e:=l to AnzEr do FZ[e]:=0; 
for v:=l to AnzVg do begin 

FAZ[v]:=FZ[i[v]]; FEZ[v]:=FAZ[v]+D[v]; 
if FEZ[v]>FZ[j[v]] then FZ[j[v]]:=FEZ[v]; 
if FEZ[v]>FE then FE:=FEZ[v]; 
end; 

{Rückwärtsrechnung} 

for e:=l to AnzEr do SZ[e]:=FE; 
for v:= AnzVg downto 1 do begin 

SEZ[v]:=SZ[j[v]]; SAZ[v]:=SEZ[v]-D[v]; 
if SAZ[v]<SZ[i[v]] then SZ[i[v]]:=SAZ[v]; 
end; 


{Tabelle der Vorgänge} 

write(' Vorgang '); 

writeln(' i j D FAZ SAZ FEZ SEZ GP'); 
for v:=l to AnzVg do begin 

write(N[v]:21,i[v]:5,j[v]:5,D[v]:5,FAZ[v]:5); 
write(SAZ[v]:5,FEZ[v]:5,SEZ[v]:5,SAZ[v]-FAZ[v]:5); 
if SAZ[v]=FAZ[v] then write(' K');writeln; 
end; 


{Tabelle der Ereignisse} 

writeln; writeln(' Ereignis FZ SZ'); 

for e:=l to AnzEr do writeln(e:17,FZ[e]:8,SZ[e]:7); readln; 

end. 
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Mit der Eingabedatei CPMPAS.DAT 

Motor ausbauen 
12 4 

Zylinderkopf zerl. 

2 3 1 

Kurbelw,Kolb.ausb. 

2 5 2 

Kolben ersetzen 

2 7 5 

Ventile einschleifen 

3 4 3 

Kopf zus.bauen 

4 8 2 

Kurbelw.schleifen 

5 6 7 

Zylinder schleifen 

5 7 6 

Pleuellager erneuern 

6 7 1 

Motorblock zus.bauen 

7 8 6 

Motor einbauen 

8 9 8 

erhalten wir die Ausgabe: 

Vorgang i 

Motor ausbauen 1 

Zylinderkopf zerl. 2 

Kurbelw,Kolben ausb. 2 

Kolben ersetzen 2 

Ventile einschleifen 3 

Kopf zus.bauen 4 

Kurbelw.schleifen 5 

Zylinder schleifen 5 

Pleuellager erneuern 6 

Motorblock zusbauen 7 

Motor einbauen 8 

FZ 
0 


8 

6 

13 

14 
20 
28 


FAZ SAZ FEZ SEZ GP 


14 15 10 


15 18 10 

18 20 10 


Ereignis 



23. Differenzengleichungen 


23.01 

a) 

(3), 

(4) b) (1), 

(3), (4) 


c) (1) 

, (3) 

23. 02 

a) 

2, 10, 34, 106, 322 


b) 1, 2, 

, 3, 

5, 8 



c) 

64, 

32, 16, 8, 4 


d) 7,5; 

15; 

22,5; 

30; 37,5 

23.03 

a) 

(1) 

konstant 


-2, 

-2, 

-2, 

-2, -2 



(2) 

fallend 


-3, 

-5, 

-11/ 

-29, -83 



(3) 

steigend 


-1, 

1 , 

7, 

25, 79 


b) 

(1) 

konstant 


0 , 

0 , 

0 , 

0 , 0 



(2) 

fallend 


1 , 

0.5 

, 0.25, 

0.125 



(3) 

steigend 


- 1 , 

-0.5 

, -0.25, 

-0.125 


c) 

(1) 

konstant 


0 , 

0 , 

0 , 

0 , 0 



(2) 

alternierend 


2, 

-2, 

2, 

-2, 2 



(3) 

alternierend 


2, 

2, 

-2, 

2, -2 


d) 

(1) 

ab a 2 fallend 


1 , 

-2, 

-1, 

-3, -4 



(2) 

ab a 5 fallend 


-2, 

1 , 

-1, 

0 , -1 

23.04 

linear: 

a, b, c, e 







konstante Koeffizienten: 


a, b, c 




23.05 

a) 

2 


b) 

5 


c) 1 



d) 

1 


e) 

2 


f) 2 


23.06 


3 

11 

27 

59 


123 

23.07 


0 

0 


3 

15 


54 

23.08 


3 

3,16 


4,2752 

7,07814 

13,0962 

23.09 


-10 

-3,1933 


-0,9244 

-o. 

1681 

0,0084 

23.10 


-o. 

3 -0,38 


-0,828 

-2, 

2968 

-6,77807 

23.11 


2 

-3,5 


15,5 

-36 

.5 

9.5 

23.12 

a n 

= C 1 

3 n + C 2 2 n 






23.13 

(-: 

L) n 3 

n (C x + C 2 n) 






23.14 

y n 

= 3 n 

cos(n arctan (4 

V2)) + C 2 

sin(n arctan (4 v'S')) j 


23.15 Z n = C x *5 n + (-1) n C 2 -4 n 



23. Differenzengleichungen 


23.16 X n+1 = 0,5 X R + 1000 x Q = 1200 

RECURRENCE1(0.5 x + 1000, n, X, 0, 1200, 10) 

[[0, 1200], [1, 1600], [2, 1800], [3, 1900], [4, 1950] 

.[5, 1975], [6, 1987.5], [7, 1993.75], [8, 1996.87], 

[9, 1998.43], [10, 1999.21]] 


23.17 X n+1 = 1,25 x n 

RECURRENCE1(1.25 X, n, X, 


X Q = 1000 
0 , 1000 , 10 ) 



23. 18 X n+1 = 0,5 X n - 100 


1800 


RECURRENCE1(0.5 X - 100, n, x, 0, 1800, 8) 







106 


23. Differenzengleichungen 


RECURRENCE1(1.05 X, n, X, 0, 250, 10) 

[[0, 250], [1, 262.5], [2, 275.625], [3, 289.406], 

[4, 303.876], [5, 319.070], [6, 335.023], [7, 351.775], 
[8, 369.363], [9, 387.832], [10, 407.223]] 


23.20 20000 (1-3/200)^60 = 8076,13 8,076 kg 

Oder LIN1_DIFFERENCE(1 - 3/200, 0, X, 0, 20000) 

20000 #e A (-0.0151136 x) 

20000 j^e A (—0.0151136 60) 

8076.13 

23.21 x_^, = 0,99 x_ = 100 nach 11 Tagen 

n+1 n 0 

LIN1_DIFFERENCE( 0.99, 0, X, 0, 100) 

100 #e A (-0.0100503 x) 

SOLVE(100 #e*(-0.0100503 x) = 90, x) 

[X = 10.4833] 

23.22 x n+1 = 1,01 x n x Q = 24000 / cm 3 nach 70 Minuten 

LIN1_DIFFERENCE(1.01, 0, X, 0, 1) 

#e A (0.00995033 x) 

SOLVE(100 #e A (0.00995033 x) = 2, X) 

[x = 69.6607] 

23.23 a) x n+1 = * n + n + 1 x Q = 2 

RECURRENCE1(x + n + 1, n, x, 0, 2, 10) 








23. Differenzengleichungen 
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RECURRENCE1( X + 2 n + 3, n, X, 0, 2, 6) 



23.24 x n+1 = x n + 500 2 n X Q = 1000 

RECURRENCE1(x + 500 2 A (-n), n, 0, 1000, 10) 

[[0, 1000], [1, 1500], [2, 1750], [3, 1875], [4, 1937.5], 
[5, 1968.75], [6, 1984.37], [7, 1992.18], [8, 1996.09], 

[9, 1998.04], [10, 1999.02]] 

23.25 x n+1 = 1,2 X n + 20 x Q = 200 

RECURRENCE1(1.2 X + 20, n, X, 0, 200, 10) 

[[0, 200], [1, 260], [2, 332], [3, 418.4], [4, 522.08], 
[5,646.496], [6, 795.795], [7, 974.954), [8, 1189.94], 

[9, 1447.93], [10, 1757.52]] 

23.26 x Q = öS 21,58 x 1 = öS 21,11 x 2 = öS 21,76 

x 3 = öS 20,86 x 4 = öS 22,10 x 5 = öS 20,39 

SOLVE([21.11 = 21.58 a + b, 21.76 = 21.11 a + b], [a, b]) 
[[a = -1.38297, b = 50.9546]] 

RECURRENCE1(-1.38297 X + 50.9546, n, X, 0, 21.58, 5) 

[[0, 21.58],[1, 21.1101], [2, 21.7599], [3, 20.8612], 

[4, 22.1041], [5, 20.3852]] 

23.27 Der Preis pendelt sich bei öS 15,14 ein. 

SOLVE([22 = 6 a + b, 10 = 22 a + b], [a, b]) 

[[a = -0.75, b = 26.5]] 

RECURRENCE1( -0.75 X + 26.5, n, X, 0, 6, 60) 

[[0, 6], [1, 22], [2, 10], [3, 19], [4, 12.25], ... 

[57, 15.1428], [58, 15.1428], [59, 15.1428], [60, 15.1428]] 
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23. Differenzengleichungen 



23.28 X n+1 = 1,2 X n - 40 l.l” X Q = 250 

RECURRENCE1(1.2 x - 40 (l.l)*n, n, X, 0, 250, 10) 

[[0, 250], [1/ 260], (2, 268], [3, 273.2], [4, 274.6], 

[5,270.956], [6, 260.726], [7, 242.009], [8, 212.462], 

[9, 169.211], [10, 108.736]] 

23.29 X R+1 = 1,085 X n - 80000 l,03 n X Q = 1200000 

RECURRENCE1(1.085 X - 80000 1.03 A n, n, X, 0, 1200000, 34) 
[[0, 1.2'10*6], [1, 1.222*10*6], [2, 1.24347•10*6], ... 

... ,[32, 2.81948-10*5],[33, 99907.8],[34, -1.03786 10*5] 

23.30 x n+1 = x n + 1000 (0.95) n - x n (25 - n)/100 

RECURRENCE1(x + 1000 (0.95)*n - x(25 - n)/100, 
n, x, 0, 5000, 20) 

5000 4750 4560 4414 4300 4212 4143 

4091 4053 4027 4013 4010 4017 4035 

4065 4105 4158 4224 4304 4400 4513 








24. Mathematik für Wirtschaftsingenieure 

24.01 

1: x := [[0, 1, 3], [2, 0, 8], [4, 10, 0]] 

2: v := [10, 20, 30] 

3: a := VECTOR(VECTOR(ELEMENT (X, i, j) / ELEMENT(W, j), j, 1, 3), 

i, 1, 3) 

4: [[0, 1/20, 1/10], [1/5, 0, 4/15], [2/5, 1/2, 0]] 

5: al := IDENTITY_MATRIX(3) - a 

6: [[1, -1/20, -1/10], [-1/5, 1, -4/15], [-2/5, -1/2, 1]] 

7: v := alW[ 10, 30, 50] [3.5, 14.6666, 31] 

9: inval := al A (-l) 

10: vneu := invalW[10, 5, 20] [14.2678, 16.9717, 34.1930] 


DERIVE - Version 3.05-G 


DERIVE verdankt seine große Verbreitung der Tatsache, daß die ersten 
Versionen auch auf einfachen Rechnern liefen. Durch die rasante 
Entwicklung der PC wurden diese immer billiger und leistungsfähiger. 
286-er und 386-er PC sind nur noch als Altgeräte zu erhalten. Nun 
verliert der 486-er seine Spitzenstellung, der Pentium ist Standard. 

Alle CA-Systeme entwickeln sich rasch weiter, die Sicherheit und die 
Bedienungsfreundlichkeit werden immer größer. 

DERIVE hat mit seinen letzten Varianten einen großen Sprung bezüglich 
der Grafik gemacht. 

Nun sind wir bei der Version 3.05 - G angelangt, die für 
deutschsprachige Benutzer gedacht ist. 

Wir bringen auf den letzten zwei Seiten eine Entsprechungsliste der 
englischen und der deutschen Menüs. 
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DERIVE Version 3 Hagenberg, im April 1995 

_ Entsprechungsliste englische Menüs <-» deutsche Menüs _ 


Author.Schreibe 

Build .Baue 

Calculus.Analysis 

Differentiate.Differenziere 

Int eg rate.Integriere 

Limit.Grenzwert 

Product.Produkt 

Sum.Summe 

Taylor.Taylor 

Vector.Vektor 

Deciare.Def 

Function.Funktion 

Variable.Variable 

Value.Value 

Integer.Integer 

Real.Real 

Complex.Complex 

Nonscalar..Nonscalar 

Matrix.Matrix 

vectoR.veKtor 

Exoand.Mult 

Factor.faKt 

Help ..Hilfe 

Editing.Eingabe 

Functions.Funktionen 

Algebra......Algebra 

2D-plot.2D-Graphik 

3D-plot.3D-Graphik 

Utility.Zusatzdateien 

State...System 

Return.Resume 

Jump .geheZu 

soLve.Löse 

Manage...zusaTz 

Annotate.Kommentar 

Branch.Wurzel 

Exponential.Exponentialumformg. en 

Logarithm. 

Logarithm. Umformungen 

Ordering.Variablenordnung 

Renumber.Numeriere_neu 

Substitute.Substituiere 

Trigonometry.Trig.umformungen 

Options.Einstellungen 

Color.. ....Farbe 

Menu.Menü 

Work...Fenster 

Display...Schirm 


Execute.Dos 

Input.Eingabe 

Mute.Warnton 

Notation.Notation 

Output.Ausgabe 

Precision.Genauigkeit 

Radix.Basis 

Plot .Graphik 

Quit .eNde 

Remove.lösche 

Simplify.Vereinfache 

Transfer.Überträge 

Load.Laden 

Derive.Derive 

State.System 

daTa...dAten 

Utility.Zusatzdatei 

Save.Speichern 

Derive.Derive 

Basic.Basic 

C.C 

Fortran.Fortran 

Pascal.Pascal 

Options.Einstellungen 

State.System 

Merge.Anhängen 

Clear.Neubeginn 

Demo.deMo 

Print.Druck 

Expressions.Ausdrücke 

Screen.Schirm 

Layout.Layout 

Options.Einstellungen 

Unremove.Rück 

moVe.schiebe 

Window..Fenster 

Close.Close 

Designate.Designate 

Flip...:.Flip 

Goto.Goto 

Next.Next 

Open.Open 

Previous.Previous 

Split.Split 

Horizontal..Horizontal 

Vertical.Vertical 

approX.approX 
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DERIVE Version 3 

Entsprechungsliste englische Menüs <-> deutsche Menüs (Graphikfenster) 


2D-Graphikfenster: 


Algebra.Algebra 

Center.zenTriere 

Delete.Lösche 

All.Alles 

Butlast.Ohnejetzten 

First.Ersten 

Last.Letzten 

Help .Hilfe 

Move .Kreuzkoordinaten 

Options.Einstellungen 

Accuracy.Zeichengenauigkeit 

Color.Farbe 

Auto.farbWechsel 

Plot.Graphik 

Menu ..Menü 

Work.Fenster 

Display.biLdschirm 

Execute.Dos 

Mute...Warnton 

Notation.Notation 

Precision.Rechengenauigkeit 

Radix.......Basis 

State...Modus 

Plot .Zeichne 

Quit ...eNde 

Range.Bildausschnitt 

Scale.Maßstab 

T ransfer.Übertrage 

Load.Laden 

State.System 

Save.Speichern 

State.System 

Print.Druck 

Screen.Schirm 

Layout..Layout 

Options.Einstellungen 

Window.Fenster 

aXes .aChsen 

Zoom.zOome 


3D-Graphikfenster: 


Algebra...Algebra 

Center..-.zenTriere 

Eye .augPunkt 

Focal .Brennpunkt 

Grids .Gitter 

Hide .Unsichtbar 

Length.Längen 

Options.Einstellungen 

Color.Farbe 

Plot.Graphik 

Menu.Menü 

Work.Fenster 

Display.biLdschirm 

Execute.Dos 

Mute.Warnton 

Notation.Notation 

Precision.Rechengenauigkeit 

Radix.Basis 

Plot .Zeichne 

Quit ..eNde 

T ransfer....Übertrage 

Acrospin.Acrospin 

Run.. Ausführen 

Save.Speichern 

Load.Laden 

State.System 

Save.Speichern 

State.System 

Print.Druck 

Screen.Schirm 

Layout.Layout 

Options.Einstellungen 

Window..Fenster 

aXes . aChsen 

Zoom.zOome 


bk teachware Lehrmittel GmbH&CoKG, Softwarepark, A-4232 Hagenberg 
Tel. +43-7236-6065, Fax +43-7236-3769 
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